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1.  INTRODUCTION 


This  user’s  guide  details  the  FANS-3D  code  model  and  the  procedure  of  execution  of  the  model. 
This  guide  was  developed  in  support  of  Environmental  Security  Technology  Certification  Program 
(ESTEP)  Project  ER-201-031. 

Dr.  Hamn-Ching  Chen  and  his  students  and  collaborators  developed  the  FANS-3D  code  over  the 
past  25  years.  Programmers  use  this  general-purpose  computational  fluid  dynamics  (CFD)  code  for 
solving  the  Navier-Stokes  equations  governing  laminar  and  turbulent  flows  in  body-fitted  curvilinear 
grids.  The  code  employs  multi-block  overset  (chimera)  grids,  including  fully  matched,  arbitrarily 
embedded,  and/or  overlapping  grids  to  facilitate  detailed  resolution  of  unsteady  laminar  and  turbulent 
flows  around  complex  geometries  involving  arbitrary  body  motions  as  well  as  fluid-structure 
interactions.  Communication  between  grid  components  is  achieved  by  Lagrange  interpolation  at  the 
fringes.  The  code  is  fully  coupled  with  the  hole-making  and  donor-finding  algorithm,  allowing  for  the 
relative  movement  of  the  grid  blocks  at  each  time  step  for  time-domain  simulation  of  fluid- 
structure  interaction  problems,  including  violent  free  surface  motions. 

The  underlying  theory  of  the  local-analytic-based  discretization  (also  known  as  finite  analytic 
based  discretization)  is  briefly  presented  in  Section  1 . 1 .  A  complete  description  of  the  formulation, 
including  the  numerical  solution  of  well-established  two-dimensional  and  three-dimensional 
benchmarks,  is  documented  in  Pontaza,  Chen,  and  Reddy  (2005).  Additional  published  work  on  the 
theory  of  the  discretization  method  is  from  Chen  and  Chen  (1984),  Chen,  Patel,  and  Ju  (1990),  and 
Chen,  Bravo,  Chen  and  Xu  (1995). 

1.1  THEORY  AND  NUMERICAL  ALGORITHM  OF  FANS  CODE 

The  authors  mentioned  above  developed  the  finite  analytic  method  for  accurate  numerical 
simulation  of  the  time-dependent  incompressible  Navier-Stokes  equations.  To  briefly  describe  the 
formulation,  consider  a  two-dimensional  domain  partitioned  into  equal-sized  non-overlapping 
elements,  De.  We  linearized  the  Navier-Stokes  equations  in  each  element  and  write: 

h 

1  ,  ~  -  d  U 

(t/0.V)£/-— VU  =  F-_ -V  p  =  L(UfP)  in De,  (1) 

Re  at 

where  3h/<9t  is  a  discrete  representation  of  the  temporal  operator  (e.g.,  a  backward  Euler 
representation)  and  V,,  is  a  discrete  gradient  operator  in  space.  Momentarily  treating  L(U,P)  as 
known  and  constant  over  the  element,  we  see  that  the  linearized  momentum  equations  are  non- 
homogeneous  advecti on-diffusion  equations. 

Treating  each  of  the  momentum  equations  as  a  transport  equation  for  the  associated  velocity 
component,  we  use  the  natural  solution  of  the  linearized  equation  as  boundary  conditions  along  the 
edges  of  the  square  element  and  solve  the  associated  equations  by  the  method  of  separation  of 
variables  to  obtain  local  analytic  interpolants  in  terms  of  unknown  neighboring  nodal  values  of  the 
velocity  components.  The  interpolant  may  be  written  as 

—  ^ '  n  _«  /  _ 

0  =  aU  +a  L(U,P)  in  Cf.  (2) 

n= 1 

The  local  analytic  interpolants  {aj*  ,  ccf  are  functions  of  the  local  velocity  field  and  respond 

analytically  to  local  flow  conditions.  In  addition,  the  interpolants  satisfy  zeroth  and  first-order 
consistency  requirements,  and  are  always  positive.  These  properties  ensure  that  spurious  energy 
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modes  are  non-existent  in  the  scheme,  and  render  it  stable  at  high  Reynolds  numbers.  Plots  of  one  of 
the  coefficients  for  different  flow  conditions  in  a  single  element  are  shown  in  Figure  1.  A  more 
detailed  description  of  the  finite  analytic  functions  is  given  in  Pontaza,  Chen,  and  Reddy  (2005)  and 
Chen  and  Chen  (1984). 


Figure  1 .  Finite  analytic  function  associated  with  a  node  placed  a 
(x,y)  =  (-1,0),  for  an  element  in  (x,y)  e  [-1,1]*  [-1,1],  for  different  flow 
conditions. 


The  interpolants  satisfy  (locally)  the  linearized  momentum  equations  and  a  collocation  scheme  is 
adopted  to  form  the  discrete  equations.  In  other  words,  the  local  analytic  functions  are  only  evaluated 
at  the  center  of  the  element  to  yield  coefficients  that  make  up  the  stencil  relating  the  center  value  to 
its  neighbors. 

If  the  pressure  field  is  known  a  priori,  the  pressure  gradient  may  be  evaluated  and  a  set  of  discrete 
equations  for  each  interior  nodded  can  be  written  using  Equation  (2).  These  equations  can  be 
assembled  to  yield  a  banded,  unsymmetrical,  definite  matrix  system.  When  augmented  with  suitable 
boundary  conditions,  the  system  can  be  solved  (in  an  iterative  manner  for  the  linearization)  to  yield 
the  nodal  velocity  values  in  a  time-marching  procedure. 

In  general,  however,  the  pressure  field  is  not  known  a  priori  and  must  be  computed  such  that  the 
velocity  field  is  divergence-free.  This  task  is  achieved  by  projecting  the  velocity  field  onto  a 
divergence-free  space  through  a  discrete  Poisson  equation  for  the  pressure.  The  discrete 
representation  of  the  divergence  operator  is  constructed  such  that  a  strong  velocity-pressure  coupling 
is  achieved,  effectively  avoiding  spurious  pressure  solutions  for  the  co-located  node  arrangement, 
where  nodal  degrees  of  freedom  for  velocities  and  pressure  share  the  same  locations.  The  projection 
is  directly  applied  to  boundaries  as  well,  so  that  no  artificial  boundary  conditions  for  the  pressure  are 
necessary.  Thus,  pressure  is  consistently  computed  at  the  boundaries. 
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The  momentum  and  discrete  pressure  Poisson  equation  are  solved  sequentially  in  an  iterative 
manner.  Pontaza,  Chen,  and  Reddy  (2005)  showed  the  method  to  be  second-order  accurate  in 
velocities  and  pressure.  Convergence  properties  of  the  method  are  illustrated  in  Figure  2.  When 
Equation  (2)  is  used  as  an  interpolant,  the  interpolation  is  fourth-order  accurate,  as  shown  in  Figure 
2(a).  When  Equation  (2)  is  used  as  a  collocation  discretization  procedure,  the  error  decays  at  a 
second-order  rate,  as  shown  in  Figures  2(b)  and  2(c)  for  linear  and  nonlinear  equations.  Figure  2(d) 
shows  second-order  accuracy  in  velocities  and  pressures,  indicating  good  velocity  pressure  coupling 
by  implementing  the  segregated  solution  approach. 

In  practical  implementations,  we  seldom  encounter  square  domains.  The  general  procedure 
consists  of  constructing  the  local  analytic  interpolants  in  a  mapped  space.  Using  this  approach,  we 
can  handle  skewed  or  curvilinear  elements  with  a  unified  approach.  The  method  has  proven  robust  in 
the  presence  of  severe  mesh  skews  and  high  aspect  ratio  cells  (Pontaza,  Chen,  and  Reddy,  2005). 


Figure  2.  Convergence  curves  for  verification  studies  of  the  finite  analytic  functions 
as  (a)  interpolants,  and  as  (b),  (c),  (d)  a  collocation  discretization  procedure. 

For  time-accurate  solutions,  the  time  derivative  is  represented  here  by  second-order  accurate 
truncated  expansions  in  time  domain.  Specifically,  the  time  integration  scheme  corresponds  to  the 
generalized  a-method  family  of  time  integrators.  The  family  is  generated  by  varying  a  single  free- 
integrator  parameter,  p,  for  high-frequency  damping.  Unresolved  high  frequencies  (due  to  the  choice 
of  the  time  step  size)  are  damped  out  according  to  the  value  of  p.  The  choice  p  =  1.0,  corresponds  to 
the  trapezoidal  rule,  which  is  known  to  have  no  damping  for  high-frequency  modes  that  may  excite 
odd-even  mode  oscillations.  High-frequency  damping  is  allowed  by  decreasing  the  value  of  p. 
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Additional  documentation  on  this  particular  family  of  time  integrators  is  given  by  Chung  and  Hulbert 
(1993)  and  by  Dettmer  and  Peric  (2003).  The  discrete  pressure  gradient  operator  is  represented  using 
standard  second-order  accurate  finite  differences  in  each  spatial  direction. 

Extension  to  the  three-dimensional  case  is  straightforward  and  achieved  by  superimposing  two- 
dimensional  local  analytic  solutions,  such  that  the  three-dimensional  equations  are  satisfied  locally. 
Details  of  the  derivation  were  first  presented  by  Chen,  Patel  and  Ju  (1990)  and  Chen,  Bravo,  Chen, 
and  Xu  (1995),  and  are  also  outlined  by  Pontaza,  Chen,  and  Reddy  (2005).  The  resulting  stencil 
relates  one  nodal  unknown  to  its  19  neighbors,  and  is  thus  a  19-point  finite  analytic  stencil. 

For  turbulent  flows  modeled  through  the  numerical  solution  of  the  Reynolds-averaged  Navier- 
Stokes  (RANS)  equations,  see  Equations  (5)  through  (11)  provided  by  Wang  and  Chen  (2016)  in 
the  technical  report  published  by  Space  and  Naval  Warfare  Systems  Center  (SSC  Pacific).  In  the 
two-layer  k-e  model,  the  k-e  model  is  patched  together  with  a  k-l  model  used  in  the  near-wall  region. 
Thus,  the  near-wall  region  is  computed  directly  and  adequate  grid  resolution  must  be  used  there. 
Additional  details  can  be  found  in  Chen  and  Patel  (1988). 

The  discretization  procedure  for  the  turbulent  transport  equations  is  exactly  the  same  used  for  the 
momentum  equations  described  earlier,  as  these  equations  can  always  be  written  in  the  standard  form 
given  by  Equation  (13).  This  standard  is  certainly  a  major  advantage  of  the  formulation,  as  no  special 
treatment  is  needed  for  the  turbulence  transport  equations. 
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2.  FANS-3D  SOFTWARE  DOCUMENTATION  AND  EXECUTION 


In  this  current  study,  the  FANS-3D  code  was  employed  for  the  propeller  wash  simulations  of  both 
DDG  5 1  ship  and  tugboat  cases  as  described  in  previous  sections.  The  computer  code  executables, 
numerical  grids,  input  files,  simulation  results,  and  animation  movies  for  all  seven  propeller  wash 
scenarios  were  delivered  to  Dr.  Pei-Fang  Wang  from  SSC  Pacific.  The  deliverables  are  organized  in 
seven  tar  (tape  archive)  files  as  follows: 

1.  ddg51_5kt_33ft.tar.gz:  DDG  51  ship  at  5  kts  and  33  ft  water  depth 

2.  ddg51_5kt_38ft.tar.gz:  DDG  51  ship  at  5  kts  and  38  ft  water  depth 

3.  ddg51_10kt_33ft.tar.gz:  DDG  51  ship  at  10  kts  and  33  ft  water  depth 

4.  ddg51_10kt_38ft.tar.gz:  DDG  51  ship  at  10  kts  and  33  ft  water  depth 

5.  tugboat_casel.tar.gz:  tugboat  scenario  1  with  propeller  blowing  to  open  water 

6.  tugboat_case2.tar.gz:  tugboat  scenario  2  with  propeller  blowing  to  pier  wall 

7.  tugboat_case3.tar.gz:  tugboat  scenario  3  with  propeller  blowing  parallel  to  pier  wall 

Each  folder  contains  the  following  set  of  files  that  must  be  written  by  the  users: 

1.  gridgen0.dat  (or  plot3d0.dat),  this  file  contains  the  multi-block  numerical  grids  in 
either  GRIDGEN  or  PLOT3D  format.  The  file  format  is  given  later. 

2.  inputblk.dat,  this  file  assigns  a  name  to  each  of  the  computational  grid  blocks  and 
contains  information  regarding  their  size  (both  active  and  phantom  grids  are  listed). 

3.  inputmpd .  dat,  this  file  contains  the  multi-processor  distribution  infomiation. 

4.  input .  dat,  this  is  the  control  program  file,  where  the  user  may  specify,  for  example,  the 
Reynolds  number,  the  time  step  size,  relaxation  factors,  etc. 

5.  overset .  in,  this  is  the  control  file  for  the  hole-cutting  and  donor-searching  program 

6.  *  .  bcs,  files  containing  the  boundary  condition  input  for  each  block  in  each  process,  a  total 
of  “number  of  processes”  files  must  be  present. 

In  all  FANS-3D  simulations,  it  is  necessary  to  construct  first  the  numerical  grid  for  each  test  case. 
The  name  of  the  grid  file  is  specified  in  inputblk .  dat.  The  grid  file  may  be  written  in  either 
GRIDGEN  or  PLOT3D  format  as  follows: 

(A)  GRIDGEN  format  (i format  =1) 

!  read  the  volume  grid  from  gridgen0.dat  file  (specified  in  inputblk.dat) 

!  each  block  has  size  nxi_GL,  net_GL,  nzt_GL 

do  nbk_GL=l , nblocks_GL  +  nphantoms_GL 
ijkst  GL=i j kpos_GL (nbk_GL)  +  1 

i j  knd_GL=i j  kpos^GL (nbk  GL) +nxi  GL (nbk_GL) *net_GL (nbk_GL) *nzt_GL (nbk_GL) 
read (10,*)  nbk  dum, nxi  GL (nbk_GL) , net_GL (nbk  GL),nzt  GL(nbk,GL) 
read(10,*)  (xref_GL ( i j k_GL) , i j k_GL=i j kst_GL, i j knd_GL) ,  & 

(yref  GL ( i j k_GL) , i j k_GL=i j kst_GL, i j knd_GL) ,  & 

( zref_GL ( i j  k  GL) , i j  k_GL=i j  kst_GL, i j  knd_GL) 

end  do 

(B)  PLOT3D  format  (i format  =  2) 

!  read  the  volume  grid  from  plot3d0.dat  file  (specified  in  inputblk.dat) 

!  each  block  has  size  nxi  GL,  net  GL,  nzt  GL 
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read ( 10 , * )  ndum 

do  nbk_GL=l , nblocks_GL  +  nphantoms_GL 

read (10,*)  nxi_GL(nbk  GL) , net_GL (nbk_GL) , nzt  GL(nbk_GL) 
end  do 

do  nbk_GL=l , nblocks  GL  +  nphantoms_GL 
ijkst  GL=i j kpos_GL (nbk_GL)  +  1 

i j  knd_GL=i j  kpos_GL (nbk_GL) +nxi  GL (nbk_GL) *net_GL (nbk_GL) *nzt_GL (nbk_GL) 
read(10,*)  (xref_GL ( i j k_GL) , i j k_GL=i j kst_GL, i j knd_GL) ,  & 

(yref  GL ( i j k_GL) , i j k_GL=i j kst_GL, i j knd_GL)  ,  & 

( zref_GL ( i j  k  GL) , i j  k_GL=i j  kst_GL, i j  knd_GL) 

end  do 

As  the  simulation  progresses  and  the  grids  move  and  rotate  with  respect  to  one  another,  the  grid 
motions  (e.g.,  ship  motion  and  propeller  rotation)  are  updated  based  on  the  reference  configuration  in 
gridgenO  .  dat  (or  plot3d0  .  dat).  In  the  above  pseudo-code  statements  nblocks_GL  and 
nphantoms_GL  are  the  number  of  active  (computational)  blocks  and  the  number  of  phantom 
blocks,  respectively;  which  were  already  read  from  inputblk .  dat.  More  details  of  the  input  files 
and  their  contents  will  be  given  in  the  following  sections,  in  the  context  of  the  example  problems. 
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3.  FANS-3D  CODE  PARALLELIZATION 


The  FANS-3D  code  is  a  general-purpose  CFD  code  allowing  for  the  numerical  solution  of  the 
Navier-Stokes  equations  governing  incompressible  flow  in  body- fitted  grids.  The  code  allows  for 
multi-block  overset  (chimera)  grids,  which  can  be  fully  matched,  arbitrarily  embedded,  and/or 
overlapping  with  each  other.  Communication  between  grid  components  is  achieved  by  Lagrange 
interpolation  at  the  fringes.  The  code  is  fully  coupled  with  the  hole-making  and  donor-finding 
algorithm,  allowing  for  the  relative  movement  of  the  grid  blocks  at  each  time  step  for  time-domain 
simulation  of  fluid-structure  interaction  problems  including  violent  free  surface  motions. 

The  FANS-3D  code  is  written  in  Fortran  90/95  standard  with  dynamic  memory  allocation  and  is 
fully  parallelized  using  Message-Passing-Interface  (MPI)  bindings.  It  employs  a  general  data 
management  strategy  that  allows  single  or  arbitrarily  large  groups  of  consecutive  or  non-consecutive 
blocks  to  be  assigned  to  different  processors.  This  strategy  enables  us  to  achieve  optimal  load 
balancing  when  dealing  with  multi-block  structured  grids  with  vastly  different  dimensions  among 
different  grid  blocks  as  shown  below. 

Given  a  multiple  block  structured  grid  with  N  blocks  of  different  sizes,  we  would  like  to  distribute 
the  workload  amongst  P  processes.  For  example,  consider  the  case  N  =  7,  as  shown  in  Figure  3. 


Figure  3.  Multiple  block  structured  grid  showing  N  =  7  blocks,  which 
are  to  be  distributed  among  P  <  7  processes. 

The  minimum  number  of  processes  allowed  in  the  parallelized  code  is  P  =  2,  and  the  maximum  for 
this  case  would  be  P  =  7,  which  would  imply  that  each  block  is  assigned  to  a  single  process.  Having 
observed  the  above  constraints,  the  code  allows  the  user  to  distribute  the  load  in  any  other  manner. 
Below  are  some  examples  (by  no  means  exhaustive)  of  valid  load  distributions,  where  we  fix  the 
number  of  available  processors: 

Example  #  1 : 

P=2 

PI: {1,2,3}  and  P2:{4,5,6,7} 

In  this  example,  process  1  is  assigned  blocks  {1,2,3}  and  process  2  is  assigned  blocks  {4, 5, 6, 7}. 


7 


Example  #2: 

P=3 

PI : { 1,4},  P2: {2,5},  and  P3:{3,6,7} 


In  this  example,  non-consecutive  numbered  blocks  are  assigned  to  different  processors.  This  is 
particularly  advantageous,  as  the  user  need  not  order  the  blocks  in  any  particular  manner  during  and 
after  the  grid  generation  process. 

The  load  distributions  should  be  such  that  the  load  is  almost  the  same  amongst  all  processes.  This 
is  not  a  requirement  in  the  code,  but  is  recommended  to  make  efficient  use  of  the  computational 
resources. 

The  information  on  load  distribution  is  read  in  through  the  file  inputmpd.dat,  and  is  as  follows 
for  example  #1  and  #2,  respectively. 


Example  #1: 

3  4  %  blocks  per  process  for  each  process 

12  3  %  global  block  numbers  for  each  process 

45  67 

Example  #2: 

2  2  3  %  blocks  per  process  for  each  process 

1  4  %  global  block  numbers  for  each  process 

2  5 

3  67 

The  above  input  is  all  that  is  needed  by  the  code  for  it  to  understand  and  schedule  the  loads  among 
the  different  processes.  In  addition,  each  processes  expects  one  boundary  condition  file,  containing 
boundary  condition  information  for  all  the  blocks  it  was  assigned.  The  fomiat  of  the  boundary 
condition  file  is  discussed  in  Appendix  A. 
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4.  COMPUTER  PLATFORMS,  COMPILATION,  AND  EXECUTION 

The  FANS-3D  code  was  tested  on  platforms  with  Linux®  as  the  operating  system,  with  Intel' 
Fortran  90/95  compilers  and  Message  Passing  Interface  Chameleon  (MPICH)  implementations. 
Specifically,  in  the  Dell'  clusters  at  Texas  A&M  Civil  Engineering  Department,  IBM®  clusters  at 
Texas  A&M  Supercomputing  Facility,  Linux"  clusters  at  U.S.  Army  Research  Laboratory  (ARL) 
High  Performance  Computing  cluster,  and  the  Cray  XE6111  (Copper)  cluster  at  Department  of 
Defense  High  Performance  Computing  Modernization  Program  (DoD  HPCMP).  The  FANS  code  and 
executable  can  be  installed  on  a  wide  variety  of  Unix  and  Linux  "  clusters  with  Message-Passing- 
Interface  (MPI)  libraries  for  parallel  computations  using  multiple  processors.  For  simplicity,  we  will 
summarize  only  the  procedures  to  compile  and  execute  the  code  on  the  Copper  cluster  at  DoD  Open 
Research  Systems  in  the  following  sections. 

The  FANS-3D  code  consists  of  18  Fortran  90  files,  each  with  a  specific  function.  A  list  of  the  files 
accompanied  with  a  brief  description  is  as  follows: 

•  main .  f  90  is  the  master  control  file  from  which  all  other  subroutines  are  called.  The  program 
follows  a  modular-style  programming  by  making  use  of  Fortran  90  modules,  which  are  invoked 
and  used  in  this  file. 

•  global .  f  90  is  where  all  global  variables  are  defined. 

•  s  f  1  ow .  f  9  0  defines  flow  parameter  variables  such  as  the  turbulence  model  coefficients. 

•  s input .  f  90  reads-in  all  the  program  control  inputs,  allocates  memory,  and  distributes  the  load 
among  processors. 

•  geocoeff.f90  computes  and  stores  the  geometric  coefficients  associated  with  a  well-defined 
transfomiation. 

•  f  acoef  f .  f  90  computes  the  19-point  stencil  finite-analytic  coefficients. 

•  moment .  f  9  0  solves  the  moment  equations  for  the  velocity  components. 

•  press  ub.f90  computes  the  19-point  stencil  for  the  pressure  Poisson  equation,  assembles  and 
solves  the  associated  system  of  equations. 

•  turbsub .  f  90  solves  the  turbulence  model  equations. 

•  sources,  f  90  computes  the  source  functions  for  the  governing  equations. 

•  boundary .  f  90  computes  and  assigns  boundary  conditions 

•  snorms  .  f  90  computes  various  metrics,  such  as  residual  norms,  outer  iteration  norms,  time 
stepping  norms,  to  establish  convergence  of  the  iterative  solution  procedure  and  time  marching 
procedure. 

•  gmotions  .  f  90  grid  motions  file  to  control  and  impose  how  the  grids  move  relative  to  each 
other  and  compute  the  grid  velocities. 

•  datamgmt .  f  90  contains  the  subroutines  for  the  multi-block  data  management. 

•  graphics  .  f  90  generates  output  files  for  visualization. 

•  sclean.f90  deallocates  memory. 

•  overset .  f  90  grid  interpolation  program  for  overset  grids. 

•  dwssub .  f  90  computes  directional  short-crested  waves. 
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The  code  is  to  be  compiled  by  linking  the  Fortran  90/95  compiler  with  a  MPI  library  or  by  using  a 
Fortran  90/95  MPI  wrapper  (e.g.,  mpiifort ,  mpif  90  or  ftn).  When  using  MPI  as  a  library,  the 
following  is  used  to  compile  the  code  on  DoD  HPCMP  Copper  cluster: 

prompt%>  module  swap  PrgEnv-pgi  PrgEnv-intel 

prompt%>  ftn  -openmp  -02  -o  fans3d.exe  {list  of  Fortran  files} 

The  code  is  simply  run  by  typing  the  following  at  the  prompt  or  giving  the  following  command  in 
the  batch-job  file  (e.g.,  for  PBS  or  LSF  queue  managers): 

prompt%>  aprun  -n  {number  of  processors}  ./fans3d.exe  >  fans3d.out 

Note  that  the  simulation  results  for  all  seven  propeller  wash  scenarios  described  earlier  can  be 
reproduced  by  uploading  the  corresponding  tar  files  to  DoD  HPCMP  Copper  cluster  and  executing 
the  following  four  commands  (using  tugboat_case3  .  tar .  gz  as  an  example): 

1 .0  Unzip  *  .  tar .  gz  file.  The  code  executable  and  input  data  files  will  be  saved  in  a  newly 
created  folder  tugboat_case3 

Prompt%>  tar  xzf  tugboat  case3.tar.gz 

2.0  Change  programming  environment  from  the  default  'pgi'  to  'intef  Fortran 

Prompt%>  module  swap  PrgEnv-pgi  PrgEnv-intel 
3.0  Switch  to  working  directory 

Prompt%>  cd  tugboat_case3 

4 . 0  Submit  job  to  the  batch  queue  (with  appropriate  project  number  in  the  job  control 

file) 


Prompt%>  qsub  submit  pbs 


5.  FANS-3D  DATA  EXPORT 


On  FANS-3D  output,  the  following  files  are  written  out  to  visualize  the  solution  using  the 
commercial  flow  visualization  software  such  as  FieldView,  TecplotR,  or  MATLAB®: 

1 .  force  .  dat,  x,  y,  and  z  forces  exerted  on  the  propeller  blades,  ship  hull  surface,  and/or  other 
solid  surfaces. 

2.  motion .  dat,  time  history  of  six-degree-of-freedom  ship  motion 

3.  overset .  out,  output  file  containing  grid  interpolation  infonnation. 

4.  f  ans3d .  out,  output  file  for  monitoring  of  convergence  history. 

5.  restart_xyz  .  dat,  instantaneous  grid  restart  file  for  continuation  run 

6.  restart_q { number }  .  dat,  instantaneous  flow  field  restart  file  for  continuation  run. 

7.  movie_x  { number }  .  dat,  three-dimensional  output  to  visualize  the  entire  grid  at  time  step 
{ number } . 

8.  movie_q  (number}  .  dat,  three-dimensional  output  to  visualize  instantaneous  velocity  and 
pressure  fields  at  time  step  { number } . 


The  force .  dat  contains  ASCII  data  files  in  column  format.  It  can  be  read  directly  into  Tecplot®1 
or  MATLAB®  or  other  compatible  software  for  2D  line  plots  of  the  (x,  y,  z)  forces  and  moments 
(with  respect  to  the  gravity  or  center  of  rotation).  For  a  problem  involving  six  degrees-of- freedom 
(heave,  sway,  surge,  pitch,  yaw,  and  roll)  motions  under  hydrodynamic  loadings,  such  as  wave  and 
current,  the  code  will  also  output  the  motion  histories  in  motion.dat  file,  which  is  also  in  ASCII 
column  data  format. 

The  overset .  out  is  an  ASCII  file  containing  grid  interpolation  infonnation  such  as 
interpolation  stencils  and  interpolation  coefficients  for  the  multi-block  overset  grid  system.  The 
f  ans3d .  out  is  also  in  ASCII  format.  It  is  used  to  monitor  the  convergence  histories  of  all  flow 
variables.  These  files  are  useful  for  debugging  of  the  input  data  files. 

The  restart  files  restart_xyz  .  dat  and  restart_q*  .  dat  are  unformatted  files  which  are 
used  internally  by  the  FANS-3D  code  for  continuation  runs.  The  code  will  automatically  read  in  the 
restart  files  if  the  users  wish  to  continue  a  previous  simulation  for  a  longer  duration. 

The  movie_x*  .  dat  and  movie_q*  .  dat  output  files  were  written  in  standard  PLOT3D  format 
as  follows: 

!  PLOT3D  grid  output  (movie  x{ number} .dat)  for  flow  visualization 
write  (54)  nblocks^GL 
write  ( 54 ) 

((nxi_GL(nbk  GL) , net_GL (nbk  GL),nzt  GL(nbk  GL) ) ,nbk_GL=l,  nblocks_GL) 

do  nbk_GL=l , nblocks_GL 

i j  kst=i j kpos_GL (nbk_GL) +1 

i j  knd=i j  kpos_GL (nbk_GL) +nxi_GL (nbk_GL) *net_GL (nbk  GL) *nzt_GL (nbk  GL) 
write(54) (xp (ijk) , ijk=ijkst, ijknd) ,  & 

(yp (ijk) , ijk=ijkst, ijknd)  ,  & 

( zp (i j k) , i j k=i j kst, i j knd) ,  & 

(iblank(ijk) , i j  k=i j  kst, i j knd) 

end  do 
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!  PL0T3D  flow  output  (movie  q { number } .dat)  for  flow  visualization 
write  (55)  nblocks  GL 
write ( 55 ) 

((nxi_GL(nbk  GL) , net_GL (nbk  GL),nzt  GL(nbk  GL) ) ,nbk_GL=l, nblocks_GL) 

do  nbk_GL=l , nblocks_GL 

i j  kst=i j kpos_GL (nbk_GL) +1 

i j knd=i j kpos_GL (nbk_GL) +nxi_GL (nbk_GL) *net_GL (nbk  GL) *nzt  GL(nbk  GL) 
write (55)  alpha, fsmach, reynolds, time 
write (55) (rho (i j k) , i j k=i j kst, i j knd) ,  & 

(rho (ijk) *u (ijk) , i jk=i j kst, ij  knd) ,  & 

(rho  (ijk)  *v(ijk)  ,ijk=ijkst,ij  knd)  ,  & 

(rho (ijk) *w (i j  k) , i j k=ij  kst, i j knd) ,  & 

(pr(ijk)  ,ijk=ijkst,ij  knd) 

end  do 


The  PLOT3D  grid  output  files  (movie_x*  .  dat)  contain  the  coordinates  (x,  y,  z)  and  blanking 
information  (iblank)  for  every  grid  point  in  the  multi-block  overset  grid  system.  The  corresponding 
flow  variables,  including  density,  momentum,  and  pressure  (p,  pu,  pv,  pw,  p)  are  stored  in  PLOT3D 
output  files  (movie  q*  .  dat).  These  data  files  can  be  imported  directly  into  the  commercial 
software  FieldView  for  flow  visualization  and  saved  in  animation  video  files  (in  avi  format).  The 
movie  data  files  can  also  be  imported  into  the  commercial  Tecplot®  software  using  the  ‘PLOT3D 
Loader’  option.  Typical  results  include  the  velocity  contours,  velocity  vector  plots,  and  pressure 
contours.  Other  quantities  such  as  shear  stresses  and  vorticities  can  also  be  calculated  using  the  user- 
defined  functions  in  FieldView  and  Tecplot®.  The  users  may  consult  the  FieldView  and  Tec)PlotR 
manuals  for  additional  information  on  the  post-processing  of  PLOT3D  data. 
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6.  EXAMPLE  CASE  1 

DDG  51  SHIP  AND  P4876  PROPELLER  WASH  STUDY 


In  this  section,  we  present  an  example  test  case  for  the  DDG  5 1  propeller  wash  study.  The  problem 
demonstrates  the  many  capabilities  of  the  FANS-3D  formulation  and  implementation,  which  include: 
embedded  and  non-matching  grids,  relative  motion  between  grid  components,  load  distribution 
among  different  processes,  high  Reynolds  number  flows,  and  robustness  in  the  presence  of  high 
aspect  ratio  skewed  meshes. 

Chen  and  Wang  (2016)  show  the  computational  domain  and  multi-block  overset  grids  for  this  case 
in  Figure  2  of  the  SSC  Pacific  technical  report.  The  length  of  the  DDG  51  ship  is  142.04  m  (466  ft) 
and  the  designed  draft  is  9.4488  m.  The  diameter  of  the  twin-screw  P4876  propellers  is  5.4864  m  (18 
ft),  and  the  center  of  propeller  axis  is  located  at  5.7912  m  below  the  mean  water  level.  We 
performed  a  calculation  for  a  shallow  water  case  with  water  depth  H  =  10.0584  m  (33  ft).  Under  this 
condition,  the  underkeel  clearance  is  only  0.6096  m  (2  ft)  beneath  the  sonar  dome  and  the  minimum 
gap  between  the  propeller  tip  and  the  sea  bottom  is  1.524  m  (5  ft).  The  twin-screw  propellers  are 
rotating  at  51  rpm  when  the  ship  speed  is  10  kts. 

A  commercial  grid  generation  software  Gridgem  was  used  to  generate  the  overset  grid  system  for 
the  DDG  51  ship  and  the  five-blade  P4876  propeller.  As  noted  earlier,  the  composite  grid  consists  of 
15  computational  blocks  and  7  phantom  grid  blocks  with  2,369,549  grid  points  covering  half  of  the 
solution  domain.  There  are  five  blocks  for  five  propeller  blades,  three  blocks  for  propeller  shaft  and 
near-wake  regions,  one  block  for  the  ship,  and  six  blocks  for  the  far  field.  The  15  blocks  are  shown  in 
different  colors  in  Figure  2  of  the  SSC  Pacific  technical  report.  In  addition,  seven  phantom  grids  (not 
shown)  are  needed  to  perform  the  hole-cutting  adequately.  The  end-user  does  not  need  to  be 
concerned  with  phantom  grids,  as  they  do  not  enter  into  the  actual  computations,  and  hence  do  not 
need  to  be  listed  in  the  multi-processor  input  file  or  the  boundary  condition  input  files. 

In  this  particular  run,  the  five  propeller  blades,  the  shaft  block,  are  assigned  to  three  processes,  the 
ship  is  assigned  to  the  fourth  process,  the  propeller  near-wake  region  is  divided  into  two  blocks  and 
assigned  to  two  separate  processes,  and  the  far-field  grids  are  decomposed  into  six  blocks  and 
assigned  to  six  different  processes.  For  this  example,  the  file  inputblk .  dat  contains  the  following 
data: 

!  Geometry  input  file  (second  line,  no  more  than  40  characters) 
gridgenO . dat 

1  !  1:  Gridgen  format,  2:  Plot3d  format 

15  7  !  nblocks  +  nphantom  (including  phantom  grid) 

62  41  41 
propellerO 1 

62  41  41 
propeller02 

62  41  41 
propeller03 

62  41  41 
propeller04 

62  41  41 
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propeller05 


38  21  122 
shaf tOl 

28  32  122 
shaf t02 

65  21  122 
shaft03 

121  35  41 
shipO 1 

34  81  77 
basinO 1 

34  81  77 
basin02 

34  81  77 
basin03 

152  65  21 
oceanO 1 

77  65  42 
ocean02 

77  65  42 
ocean03 

3  41  61 
phantomO 1 

3  41  61 
phantom02 

3  41  61 
phantom03 

3  41  61 
phantom04 

3  41  61 
phantom05 

2  2  2 
phantomO  6 

2  2  2 
phantom07 

This  input  specifies  that  the  name  of  the  composite  grid  file  is  gridgenO  .  dat,  and  it  is  in 
GridgenR  format.  There  are  15  computational  blocks  and  7  phantom  blocks  (22  blocks  total).  Then, 
for  each  of  the  15  computational  blocks,  we  must  specify  their  (i,  j,  k)  sizes  and  assign  to  a  name  to 
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them  that  must  be  consistent  with  the  names  used  in  the  overset .  in  input  for  the  hole-cutting  and 
donor-search  algorithm. 

The  file  overset .  in  contains  the  input  necessary  for  the  hole-cutting  and  donor-search  program. 
The  format  of  this  file  is  not  discussed  here,  and  the  interested  reader  may  consult  the  Chimera 
Overset  Structured  Mesh-Interpolation  Code  (COSMIC)  Users’  Manual  (Chen,  2009).  The  input  file 
used  for  this  case  is  shown  in  Appendix  B. 

The  file  inputmpd.dat  contains  the  information  necessary  for  the  code  to  distribute  the  load 
among  the  different  processes,  as  described  in  the  previous  section.  For  this  particular  case,  the  file 
has  the  following  information: 

%  number  of  blocks  per  process,  for  each  process 
222111111111 

12  %  global  block  number  per  process,  for  each  process 

3  4 
5  6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Note  that  only  active  (computational)  blocks  are  listed  in  this  input,  i.e.,  phantom  blocks  do  not 
need  to  be  distributed  as  they  do  not  represent  any  computational  load.  In  this  particular  case,  we 
assign  propeller  blades  1  and  2  (propel lerO  1 ,  propeller 02)  to  first  process,  blades  3  and  4 
(propel ler03,  propeller 04)  to  the  second  process,  and  blade  5  and  the  first  shaft  block 
(propeller05,  shaftOl)  to  the  third  process.  The  remaining  nine  computational  blocks  (two 
shaft  blocks,  one  ship  block,  three  basin  blocks,  and  three  ocean  blocks)  are  assigned  to  processes 
#4-#  12  with  only  one  single  block  in  each  process. 
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The  input .  dat  file  is  the  main  control  input  file  and  is  as  follows: 

1  %  MTURB  flag  for  laminar  (0)  or  turbulent  (1)  flow 

1  %  INCOMP  flag  for  incompressible  (1)  or  compressible  (0)  flow 

0  %  IFSURF  flag  for  (1)  free  surface  flow  (0)  no  free  surface 

2.1868E7  %  RE  Reynolds  number 

0.04  %  TAU  time  step  size 

0.0  %  AMP  RHO  frequency  damping  parameter:  0.0  <=  AMP  RHO  <=  1.0 

1.0E-08  %  TOL1  L2  vel  tol  to  stop  time  stepping 

1.0E-03  %  TOL2  LI  res  tol  to  stop  outer  iterations 

1  %  ITIMEST  starting  time  step  to  compute 

12500  %  ITIMENDending  time  step  to  compute 

1  %  MAXIT  LS  max  allowable  ADI  sweeps  for  level-set  function 

3  %  MAXITER  max  allowable  outer  iterations 

2  %  MAXSWP  U  max  allowable  number  of  momentum  eqns  ADI  sweeps 

2  %  MAXSWP  PR  max  allowable  number  of  pressure  eqn  ADI/SIP  sweeps 

2  %  MAXSWP  KE  max  allowable  number  of  k-epsilon  eqns  ADI  sweeps 

6  %  MAXIT  DIVU  max  projections  of  velocity  field  onto  div-free  space 

0.60  %  RFU  relaxation  factor  for  velocities  (due  to  nonlinearity) 

0.30  %  RFP  relaxation  factor  for  pressure  (due  to  u-p  decoupling) 

0.010  %  RFKE  relaxation  factor  for  turbulent  k.e.  and  dissipation 

0.5  %  RFPHI  relaxation  factor  for  level-set  function 

0  %  ITIME  BCS  flag  to  indicate  (1)  time  dependent  bcs 

fans.grd  %  GEOFILE  geometry  input  file  (HCC:  not  used  in  this  version) 
ddg.bcs  %  BCSFILE  boundary  conditions  input  file 

1  %  IACT  PLOT  flag  to  activate  (1)  visualization  output 

100  %  ISKP  PLOT  time  intervals  for  vis  and  restart 

1  %  IACT  ANIME  flag  to  activate  (1)  animation  output 

0  %  IBGN  ANIME  time  step  number  at  which  animation  begins 

4  %  ISKP  ANIME  multiples  at  which  sol  is  written  out  for  animation 

0  %  ISOL  PR  (O)TDMA-ADI,  (l)SIP-7pt  solver  for  pressure  eqn 

-5.0  %  UMIN 

5.0  %  UMAX 

-10.  %  PMIN 

10.  %  PMAX 

0.050  %  TKEMAX 

0.02  %  TVISMAX 

-0.5  %  PHILSMIN 

0.5  %  PHILSMAX 

1.0  %  FROUDE  Froude  number  (gravity  acts  in  negative  z-direction) 

0.0020  %  EWIDE  representative  grid  size 

0.0  %  RFG  geometry  distortion  relaxation  parameter 

0.0  0.0  0.0  %  UINF,  VINF,  WINF  inflow  velocities 

-1.103143  0.  0.  %  UBODY,  VBODY,  WBODY  body  velocities  (ship  speed) 

1  %  IMOVE  (0)  fixed  grid,  (1)  moving  grid 

000  %  NBODY  (#bodies),  NFBODY  (fsurfaces) ,  MBLK  (#blocks  for  6-dof) 
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1  12  %  NPROP  ( #propellers ) ,  NMVPROP  (max  #blocks  moving  with  propeller) 

11%  IPROP  (propeller  ID),  IROT  (1:  clockwise,  -1:  counter-clockwise) 

12  6  %  IMVPROP  (no.  of  moving  propeller  blocks),  ISHAFT  (shaft  ID) 

1  2  3  4  5  6  7  16  17  18  19  20  %  moving  propeller  block  ID  (propeller  #1) 

Most  of  the  entries  of  the  input  file  are  self-descriptive,  but  we  further  elaborate  on  each  of  them  in 

the  following: 

•  MTURB  is  a  flag  to  specify  whether  to  numerically  solve  the  Navier-Stokes  equations  directly 
(laminar  flow  or  DNS)  or  to  solve  the  Reynolds-averaged  Navier-Stokes  equations  with  the  near¬ 
wall  two-layer  k-epsilon  model. 

•  INCOMP  is  a  flag  to  specify  whether  flow  is  incompressible  or  compressible. 

•  IFSURF  is  a  flag  to  specify  whether  it  is  necessary  to  update  free  surface.  For  this  case,  the  free 
surface  effect  is  ignored. 

•  RE  is  the  Reynolds  number.  For  this  case,  it  is  based  on  the  propeller  diameter  and  the  propeller 
rotating  speed. 

•  TAU  is  the  value  for  the  time  step  size.  In  this  case  the  angular  velocity  of  the  propeller  is  such 
that  one  revolution  is  completed  in  one  unit  of  time. 

•  AMP_RHO  is  the  high-frequency  damping  parameter  for  the  second-order  accurate  family  of  time 
integrators,  as  described  in  the  formulation  section. 

•  TOL1  is  tolerance  of  L2  velocity-norm  to  stop  time  stepping  when  the  steady  state  is  reached. 

•  TOL2  is  tolerance  of  LI  velocity-residuals  to  stop  outer  iterations  within  each  time  step. 

•  ITEMEST  is  the  starting  time  step  of  the  computation.  A  value  of  1  is  specified  for  new  runs.  If 
the  value  is  greater  than  1 ,  the  code  will  read-in  restart  files  from  previous  runs  and  continue  the 
computation  to  the  new  ending  time  step. 

•  itemend  is  the  ending  time  step  the  user  wishes  to  compute,  for  the  previously  specified  value 
of  the  time  step  size. 

•  MAXIT_LS  is  the  maximum  allowable  Alternating-Directional-Implicit  (ADI)  sweeps  for  the 
level-set  function. 

•  MAXITER  is  the  maximum  allowable  number  of  outer  iterations  on  a  given  time  step.  For  time 
accurate  solutions  this  value  must  be  greater  than  one,  to  allow  for  good  velocity-pressure 
coupling  and  hence  time  accuracy  of  the  flow  field. 

•  MAXSWP_U,  MAXSWP_P,  MAXSWP_KE  is  the  maximum  allowable  number  of  inner  iterations 
on  a  given  outer  iteration,  to  iteratively  solve  the  momentum,  pressure,  and  turbulence  transport 
equations,  respectively. 

•  MAX  I  T_D  I VU  is  the  maximum  allowable  number  of  projections  of  the  velocity  field  onto  a 
divergence-free  space  on  a  given  outer  iteration. 

•  RFU  is  the  relaxation  factor  for  the  velocity  field.  The  optimal  values  lie  in  the  range  [0.4,  1.0], 
although  lower  values  may  be  needed  for  complex  problems. 

•  RFP  is  the  relaxation  factor  for  the  pressure  field.  Also,  in  accordance  with  well-established 
practices,  we  find  that  optimal  values  lie  in  the  range  [0.1,  0.8],  although  higher  values  may  also 
be  used  and  lower  values  may  also  be  needed. 

•  RFKE  is  the  relaxation  factors  for  the  turbulent  transport  variables.  We  find  that  typically  optimal 
values  lie  in  the  range  [0.01,  0.5].  Although  lower  values  may  be  needed. 
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•  RFPHI  is  the  relaxation  factors  for  the  level-set  function.  We  find  that  typically  optimal  values 
lie  in  the  range  [0.2,  1.0]. 

•  BCSFILE  is  a  string  specifying  the  name  of  the  family  of  boundary  conditions  files.  The  family 
must  have  “number  of  processors”  members.  In  this  particular  case,  since  there  are  12  processes 
involved,  we  must  have  12  files  (ddgO  0  .  bcs  ~  ddgl  1 .  bcs)  ready. 

•  GEOFILE  is  a  string  specifying  the  name  of  the  grid  file.  It  is  not  needed  in  this  version  since  the 
grid  name  has  already  been  specified  in  inputmpd .  dat. 

•  IACT_PLOT ,  I  SKP_PLOT  are  control  flags  to  activate  the  output  and  to  control  how  frequently 
the  output  files  are  updated.  The  output  is  in  PLOT3D  format  for  visualization  using  commercial 
software  such  as  FieldView,  Tecplot®  or  other  compatible  flow  visualization  tools. 

•  IACT_ANIME ,  IBGN_ANIME ,  I  SKP_ANIME  are  control  flags  to  write  out  a  movie,  which  is 
to  be  processed  by  the  software  FieldView.  For  this  case,  the  movie  corresponds  to  grid 
coordinates  and  grid  blanking  values,  density,  velocity  vectors,  and  pressure  on  every  point  in  the 
flow  field. 

•  ISOL_PR  is  a  flag  for  pressure  solvers.  The  pressure  can  be  solved  using  either  tridiagonal 
matrix  algorithm  (TDMA-ADI)  or  strongly- implicit  method  (SIP-7pt). 

•  UMIN,  UMAX,  PMIN,  PMAX,  TKEMAX,  TVISMAX  are  limiters  on  the  velocity,  pressure, 
and  turbulent  transport  variables.  These  are  set  to  high  values,  and  are  just  a  safeguard  against  a 
poor  initial  guess,  which  may  cause  the  fields  to  oscillate  violently  in  the  initial  stages  of  the 
iterations. 

•  PHILSMIN,  PHILSMAX  are  limiters  for  the  level-set  function.  They  are  set  to  high  values,  and 
are  just  a  safeguard  against  a  poor  initial  guess,  which  may  cause  the  fields  to  oscillate  violently 
in  the  initial  stages  of  the  iterations. 

•  EWIDE  is  a  representative  grid  size  used  specifying  the  transitional  zone  thickness  adjacent  to 
the  air-water  interface. 

•  RFG  is  a  geometry  distortion  parameter.  The  default  value  is  1  for  orthogonal  or  nearly- 
orthogonal  grids,  but  may  be  reduced  to  improve  convergence  for  highly-skewed  grids.  The 
relaxation  parameter  does  not  affect  accuracy  for  orthogonal  grids  and  has  negligible  effects  for 
nearly  orthogonal  grids. 

•  UINF,  VINF,  WINF  are  the  values  of  (x,y,z)  components  of  the  free-stream  velocity.  For  this 
case,  the  ambient  current  velocity  is  zero. 

•  UBODY ,  VBODY,  wbody  are  the  values  of  (x,y,z)  components  of  the  body  velocity  (i.e.,  ship 
speed)  normalized  by  the  characteristic  velocity  nD,  when  n  is  the  propeller  rotating  speed  and  D 
is  the  propeller  diameter.  For  this  case,  the  ship  is  traveling  in  negative-x  direction  with  a 
normalized  speed  equals  to  the  propeller  advance  coefficient  J=V/nD. 

•  IMOVE  is  a  flag  for  grid  motion.  A  value  of  0  is  specified  for  fixed  grid  system.  The  value  is  set 
to  1  for  moving  grid  in  this  case  since  the  ship  is  moving  at  constant  forward  speed  and  the 
propeller  is  also  turning. 

•  NBODY ,  NFBODY,  MBLK  are  the  number  of  body  for  force/moment  integration,  maximum 
number  of  surfaces  for  force/integration,  and  the  maximum  number  of  blocks  with  six-degree-of- 
freedom  motions.  These  parameters  are  not  needed  for  the  propeller  wash  study  considered  in  the 
present  study. 

•  NPROP  is  the  number  of  propellers.  A  value  of  1  is  specified  for  single-screw  propeller.  The 
value  is  set  to  2  for  twin-screw  propellers.  In  this  case,  we  set  NPROP  =  1  since  the  computation 
was  performed  for  only  one-half  of  the  solution  domain.  A  value  of  2  should  be  specified  for 
fully  domain  calculations  involving  twin-screw  propellers. 
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•  NMVPROP  is  the  maximum  number  of  grid  blocks  rotating  with  any  propellers. 

•  I  PROP  is  the  propeller  ID.  For  twin-screw  propellers,  the  propeller  rotating  directions  and  the 
computational  grid  blocks  associated  with  each  propeller  can  be  defined  separately. 

•  I  ROT  is  a  flag  specifying  the  propeller  rotating  direction.  A  value  of  1  indicates  that  the  propeller 
is  rotating  in  clockwise  direction.  For  counter-rotating  twin-screw  propellers,  it  is  convenient  to 
straightforward  to  specify  I  ROT  =  -1  for  the  second  propeller  rotating  in  counter-clockwise 
direction.  Also,  it  is  convenient  to  change  the  signs  of  I  ROT  if  the  same  propellers  are  under  the 
crash-astern  condition. 

•  IMVPROP  is  the  total  number  of  grid  blocks  (including  phantom  grids)  rotating  with  a  given 
propeller.  In  this  case,  there  are  5  propeller  blade  blocks,  2  shaft/hub  blocks,  and  5  phantom 
blocks  (one  for  each  propeller  blade)  rotating  with  each  propeller.  The  identification  numbers  of 
the  rotating  grid  blocks  are  specified  in  the  next  line. 

•  I  SHAFT  is  the  block  identification  number  of  the  shaft  grid  block.  This  allows  the  users  to 
specify  the  center  or  rotation  for  each  propeller. 

Since  the  workload  is  distributed  to  12  processors  for  parallel  execution,  it  is  necessary  to  write  12 
separate  boundary  condition  files  (ddgO 0  .  bcs~ddgl  1 .  bcs)  which  are  included  in  the 
ddg51_10kt_33f  t .  tar .  gz  for  the  present  case.  The  boundary  condition  files  follow  the  format 
outlined  in  Appendix  A.  Below  we  discuss,  as  an  example,  the  boundary  conditions  specified  for  one 
of  the  propeller  blades. 


propellerO 1 

1  6  62  41  41  0 


%  global  block  #1 
%  mb, nf abcs , ni , n j , nk 


2 

113  1 

2  2  0  1 


%  No .  of  two-layer  regions 
62  1  21  1  41  %  nreg, iedy, idist, (i,j,k) 

62  20  41  1  41  %  nreg,  iedy, idist,  (i,j,k) 


1  %  No.  of  free  surface  regions 

1  3  1  62  1  41  1  41  %  nLSreg,  nLS,  (i,j,k) 


1  1 

4  4  4  4  4 

1  41  1  41 

2  1 

4  4  4  4  4 

1  41  1  41 

3  1 

9  9  9  2  1 

1  62  1  41 

4  1 

4  4  4  4  4 

1  62  1  41 

5  1 

9  9  9  2  1 

I  62  1  41 

6  1 

II  11  11  11 

1  62  1  41 

14  2  64  41 


%  Face  #1 

4  4  %  (u,  v,  w, p, k, epsilon, phiLS) 

%  ( i , j , k)  range 
%  Face  #2 

4  4 

%  Face  #3 

3  3 

%  Face  #4 

4  4 

%  Face  #5 

3  3 

%  Face  #6 
11  11  11 


%  nbk_GL, i , j , k  for  pressure  datum 
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The  first  line  lists  the  blockname  for  block  identification.  The  second  line  specifies  that  this  is 
block  #1  for  the  given  process,  it  has  six  faces  with  boundary  conditions,  and  the  (i,  j,  k)  dimensions 
of  the  block  are  62  x  4 1  x  4 1 .  The  fourth  line  indicates  that  two  regions  need  to  be  identified  to  apply 
the  near-wall  two-layer  k-epsilon  turbulence  model.  The  following  two  lines  specify  the  two-layer 
model  types  (iedy,  near-wall  or  outer),  the  identification  (idist)  of  each  wall  boundary,  and  the 
(i,  j ,  k)  range  of  the  specific  region.  The  eighth  line  indicates  that  only  one  region  exists  for  the  level- 
set  function  (phiLS)  specification.  The  free  surface  solver  option  (nLS)  and  the  (i,  j,  k)  range  for  that 
region  are  specified  in  the  next  line.  These  free  surface  boundary  conditions  are  not  used  in  this  case 
since  the  free  surface  effect  is  neglected  with  the  flag  IFSURF  =  0. 

Then,  for  each  face,  we  read  the  face  number  and  the  number  of  sections  in  the  face.  For  each 
section  on  a  given  face,  we  read  7  boundary  conditions  associated  with  each  of  the  7  field  variables: 

(u,  v,  w,  p,  k,  s,  phiLS),  and  the  surface  limits  on  that  face. 

Faces  #1  and  #2  (i-min  and  i-max,  respectively)  of  the  blade  block  receive  interpolation 
information,  and  thus  all  the  field  variables  have  boundary  condition  #4.  Face  #3  (j-min)  is  the  solid- 
surface  of  the  blade,  for  which  (u,  v,  w)  are  assigned  the  grid  velocity  due  to  the  rotation  of  the  blade, 
p  is  linearly  extrapolated,  turbulent  kinetic  energy  is  zero  on  the  wall,  and  the  Neumann  boundary 
conditions  are  used  for  turbulent  energy  dissipation  and  level-set  function.  Face  #4  (j-max)  receives 
interpolation  information  and  all  field  variables  have  boundary  condition  #4.  Face  #5  (k-min)  is  part 
of  the  shaft’s  solid-surface  and  its  boundary  conditions  are  identical  to  those  on  face  #3  for  a  solid- 
wall.  Face  #6  (k-max)  is  a  branch  cut  around  the  blade  tip  where  the  flow  variables  are  updated  by 
averaging  the  adjacent  nodal  values  on  either  side  of  the  branch  cut  plane. 

For  this  example  run,  the  flow  field  is  initialized  with  calm  water  condition  and  the  propeller  is 
allowed  to  rotate  for  100  revolutions.  The  ship  travels  at  a  constant  forward  speed  of  10  kts  and  the 
propeller  rotating  speed  is  51  rpm  (0.85  rps).  The  flow  conditions  correspond  to  an  advance 
coefficient  J=  1.1 03  and  a  Reynolds  number  of  2. 1868  x  107  based  on  the  propeller  diameter.  This 
corresponds  to  a  Reynolds  number  of  6.245  x  10s  based  on  the  ship  length  and  ship  speed.  The 
RANS  equations  are  solved  with  the  near- wall  two-layer  k-e  turbulence  model. 

As  noted  earlier,  the  PLOT3D  grid  output  files  (movie_x*  .  dat)  contain  (x,  y,  z,  iblank)  for  the 
multi-block  overset  grids,  while  the  PLOT3D  flow  output  files  (movie_q*  .  dat)  contain  the  flow 
variables  (p,  pu,  pv,  pw,  p).  These  data  files  can  be  imported  directly  into  FieldView  for  flow 
visualization  and  post-processing.  Typical  results  include  the  velocity  contours,  velocity  vector  plots, 
and  pressure  contours  as  shown  in  Figures  2  through  7.  Other  derived  quantities  such  as  shear 
stresses  and  vorticities  can  also  be  calculated  using  the  user-defined  functions  in  FieldView.  Figure  8 
shows  the  shear  stresses  on  the  sea  bed  which  can  be  readily  obtained  by  evaluating  the  velocity 
gradients  adjacent  to  the  bottom  boundary  using  the  following  formula: 

dq  A  q  ^ 

r  =  /u _ =  M _  ,  q=sju  +  v  +  w/ ,  (14) 

on  An  wall 

where  q  is  the  velocity  magnitude,  //  is  the  dynamic  viscosity  of  the  seawater,  and  An  is  the  normal 
distance  from  the  wall.  The  same  shear  stress  data  can  also  be  plotted  using  another  commercial  code 
MATLAB B’,  as  shown  in  Figure  9. 
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7.  EXAMPLE  CASE  2 

TUGBOAT  AND  DUCTED  PROPELLER  WASH  STUDY 


In  the  second  test  case,  we  consider  a  tugboat  boat  with  two  ducted  propellers  under  bollard-pull 
condition  (i.e.,  zero  forward  speed),  as  shown  in  Figure  1 1  of  the  SSC  Pacific  technical  report  and 
Figures  1  and  2  in  this  document.  The  composite  grid  was  generated  by  the  commercial  grid 
generation  software  Gridgen®.  It  consists  of  47  computational  blocks  and  9  phantom  grid  blocks, 
with  7,070,832  total  grid  points.  There  are  four  blades  for  each  propeller,  and  each  blade  is  divided 
into  two  overlapping  computational  blocks.  Each  ducted  propeller  assembly  is  surrounded  by  five 
computational  blocks  covering  the  upstream,  downstream,  inner,  and  outer  regions  between  the 
propeller  shaft  and  the  shroud.  In  addition,  two  near-wake  cylindrical  grid  blocks  (one  for  each 
propeller)  are  added  to  provide  accurate  resolution  of  the  propeller  wake  flows.  The  tugboat  is 
surrounded  by  a  single  boundary- fitted  grid  block,  and  the  far  field  is  covered  by  1 8  overlapping 
rectangular  grid  points.  A  near- wall  spacing  of  10"6  ft  was  used  near  the  sea  bottom  to  provide 
accurate  resolution  of  the  turbulent  boundary  layer  flow.  This  allows  us  to  calculate  the  shear  stresses 
on  the  seabed  directly  without  relying  on  the  wall-function  approximations. 

The  composite  grid  load  is  now  distributed  among  35  processes  and  we  consider  the  bollard-pull 
(zero  tugboat  speed)  condition  with  the  ducted  propellers  blowing  parallel  to  a  pier  wall.  The  file 
inputblk .  dat  contains  the  following  data: 

!  Geometry  input  file  (second  line,  no  more  than  40  characters) 
gridgenO . dat 

1  !  1:  Gridgen  format,  2:  Plot3d  format 

47  9  !  nblock  +  nphantomO  (including  phantomO  grid) 

62  35  42 
propellerO la 

62  35  42 
propeller02a 

62  35  42 
propeller03a 

62  35  42 
propeller04a 

29  4  5 

tipO la 

29  4  5 

tip02a 

29  4  5 

tip03a 

29  4  5 

tip04a 

21  57  122 
ductO la 

66  24  122 
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duct02a 


53  35  122 
duct03a 

53  23  122 
duct04a 

21  79  122 
duct05a 

40  41  122 
wakeOla 

62  35  42 
propellerO lb 

62  35  42 
propeller02b 

62  35  42 
propeller03b 

62  35  42 
propeller04b 

29  4  5 

tipOlb 

29  4  5 

tip02b 

29  4  5 

tip03b 

29  4  5 

tip04b 

21  57  122 
ductOlb 

66  24  122 
duct02b 

53  35  122 
duct03b 

53  23  122 
duct04b 

21  79  122 
duct05b 

40  41  122 
wakeOlb 

20  116  95 
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oceanOl 

20  116  95 
ocean02 

20  116  95 
ocean03 

20  116  95 
ocean04 

20  116  95 
ocean05 

20  116  95 
oceanO  6 

20  116  95 
ocean07 

20  116  95 
oceanO  8 

20  116  95 
oceanO  9 

20  116  95 
oceanlO 

20  116  95 
oceanl 1 

20  116  95 
oceanl2 

20  116  95 
oceanl3 

20  116  95 
oceanl4 

20  116  95 
oceanl5 

16  116  95 
oceanl 6 

151  37  33 
oceanl7 

151  37  33 
oceanl 8 

107  34  61 
bargeOl 

3  31  41 
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phantomOla 


3  31  41 
phantom02a 

3  31  41 
phantom03a 

3  31  41 
phantom04a 

3  31  41 
phantomOlb 

3  31  41 
phantom02b 

3  31  41 
phantom03b 

3  31  41 
phantom04b 

2  3  2 

phantom05 

The  file  inputmpd.dat  contains  the  information  necessary  for  the  code  to  distribute  the  load 
among  the  35  different  processes.  For  this  particular  case,  the  file  has  the  following  information: 

%  number  of  blocks  per  process,  for  each  process  (excluding  phantom  grids) 
44111111441111111111111111111111111 

1526  %  global  block  number  per  process,  for  each  process 

3  7  4  8 
9 

10 

11 

12 

13 

14 

15  19  16  20 
17  21  18  22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


24 


36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 


Once  again,  the  only  active  (computational)  blocks  are  listed  in  this  input,  i.e.,  phantom  blocks  do 
not  need  to  be  distributed  as  they  do  not  represent  any  computational  load.  In  this  particular  case, 
there  are  47  actives  blocks,  including  14  blocks  for  each  ducted  propeller  assembly,  2  blocks  for 
propeller  near- wakes,  1  block  for  the  tugboat,  and  18  blocks  for  the  far  field.  To  balance  the 
workload  for  each  process,  we  assign  two  blade  and  two  tip  blocks  to  a  single  processor.  More 
specifically,  the  first  eight  computational  blocks  (#1  though  #8)  consist  of  four  blade  and  four  tip 
grids  for  the  first  ducted  propeller  are  assigned  to  processes  #0  (master  process)  and  #1,  while  the 
other  eight  blade/tip  blocks  (#15  through  #22)  for  the  second  ducted  propeller  are  assigned  to 
processes  #8  and  #9  as  shown  in  the  inputmpd .  dat  file.  The  remaining  31  computational  blocks 
containing  the  ship,  shroud,  near- wake,  and  far- field  grids  are  assigned  to  3 1  different  processes.  As 
noted  earlier,  the  user  need  not  order  the  blocks  in  any  particular  manner  during  and  after  the  grid 
generation  process. 

The  input .  dat  file  (the  main  control  input  file)  requires  only  minor  modifications,  relative  to 
the  DDG  5 1  ship  case.  The  most  notable  difference  is  that  there  are  two  co-rotating  propellers  in  the 
present  full-domain  simulation.  For  each  ducted  propeller,  there  are  12  rotating  grid  blocks  (4  blade 
surface  blocks,  4  blade  tip  blocks,  and  4  phantom  grids).  It  is  straightforward  to  specify  the  rotating 
direction,  shaft  block  ID,  and  the  IDs  of  rotating  grid  blocks  associated  with  each  propeller  in  the 
following  input .  dat  file. 


1  %  MTURB 

1  %  INCOMP 

0  %  IFSURF 

2.6468E5  %  RE 

0.04  %  TAU 

0.0  %  AMP_RHO 

1.0E-08  %  TOL1 
1.0E-03  %  TOL2 

1  %  ITIMEST 


12500 
1 
6 
2 
2 
2 
2 
0, 

0, 

0, 

0, 


4 
2 

001 

5 


turbulence  model:  (0) laminar  (l)k-epsilon  (2)LES 
flag  for  incompressible  (1)  or  compressible  (0)  flow 
flag  for  free  surface  flow  (1)  or  no  free  surface  (0) 
Reynolds  number  (L=lft,  U=lft/s,  T=L/U=ls) 
time  step  size 

frequency  damping  parameter:  0.0  <=  AMP  RHO  <=  1.0 
L2  vel  tol  to  stop  time  stepping 
LI  res  tol  to  stop  outer  iterations 
starting  time  step  to  compute 
%  ITIMENDending  time  step  to  compute 

MAXITER  LS  max  allowable  outer  equation  for  level-set  eqn 
MAXITER  max  allowable  outer  iterations 

MAXSWP  U  max  allowable  number  of  momentum  eqns  ADI  sweeps 
MAXSWP  PR  max  allowable  number  of  pressure  eqn  ADI/SIP  sweeps 
MAXSWP  KE  max  allowable  number  of  k-epsilon  eqns  ADI  sweeps 
MAXIT  DIVU  max  projections  of  velocity  field  onto  div-free  space 
;  RFU  relaxation  factor  for  velocities  (due  to  nonlinearity) 

;  RFP  relaxation  factor  for  pressure  (due  to  u-p  decoupling) 

%  RFKE  relaxation  factor  for  turbulent  k.e.  and  dissipation 
%  RFPHI  relaxation  factor  for  level  set  funtion 
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0  %  ITIME  BCS  flag  to  indicate  (1)  time  dependent  bcs 

fans.grd  %  GEOFILE  geometry  input  file  (HCC:  not  used  in  this  version) 
prop. bcs  %  BCSFILE  boundary  conditions  input  file 

1  %  IACT  PLOT  flag  to  activate  (1)  visualization  output 

100  %  ISKP  PLOT  time  intervals  for  vis  and  restart 

1  %  IACT  ANIME  flag  to  activate  (1)  animation  output 

0  %  IBGN  ANIME  time  step  number  at  which  animation  begins 

20  %  ISKP  ANIME  multiples  at  which  sol  is  written  out  for  animation 

0  %  ISOL  PR  (O)TDMA-ADI,  (l)SIP-7pt  solver  for  pressure  eqn 

-30.  %  UMIN 

30 .  %  UMAX 

-200.  %  PMIN 

200.  %  PMAX 

1.0  %  TKEMAX 

0.01  %  TVISMAX 

-30.  %  PHILSMIN 

30.  %  PHILSMAX 

1.0  %  FROUDE  Froude  number  (gravity  acts  in  negative  z-direction) 

0.003  %  EWIDE  representative  grid  size 

0.0  %  RFG  geometry  distortion  relaxation  parameter 

0.0  0.0  0.0  %  UINF,  VINF,  WINF  inflow  velocities 

0.0  0.0  0.0  %  UBODY,  VBODY,  WBODY  body  velocities  (ship  speed) 

1  %  IMOVE  (0)  fixed  grid,  (1)  moving  grid 

000  %  NBODY  (#bodies),  NFBODY  (fsurfaces) ,  MBLK  (#blocks  for  6-dof) 

2  12  %  NPROP  ( #propellers ) ,  NMVPROP  (max  #blocks  moving  with  a  propeller) 

11  %  IPROP  (propeller  #1),  IROT  (1:  clockwise,  -1:  counter-clockwise) 

12  10  %  IMVPROP  (no.  of  moving  propeller  blocks),  ISHAFT  (shaft  ID) 

1  2  3  4  5  6  7  8  48  49  50  51  %  moving  propeller  block  ID  (propeller  #1) 

21  %  IPROP  (propeller  #2),  IROT  (1:  clockwise,  -1:  counter-clockwise) 

12  24  %  NMVPROP  (no.  of  moving  propeller  blocks),  ISHAFT  (shaft  ID)  15 

16  17  18  19  20  21  22  52  53  54  55  %  moving  propeller  ID  (propeller  #2) 

The  file  overset .  in  needs  to  be  suitably  modified  for  the  hole-cutting  and  donor-search 
algorithm,  and  is  provided  in  Appendix  C.  New  boundary  condition  files  need  to  be  created  for  the 
new  tugboat  and  ducted  propeller  geometries.  In  addition,  the  boundary  condition  files  for  the  far 
field  need  to  be  modified  slightly  to  enforce  the  no-slip  boundary  conditions  on  pier  walls.  All  35 
boundary  condition  files  are  included  in  tugboat_case3  .  tar .  gz  zipped  folder. 

For  this  example  run,  the  flow  field  is  initialized  with  a  calm  water  condition  and  the  propeller  is 
allowed  to  rotate  for  500  revolutions  under  the  bollard-pull  condition  with  zero  forward  speed.  The 
simulation  was  performed  for  12,500  time  steps  with  a  time  increment  of  0.04  T0,  where  T0  is  a 
characteristic  time  for  the  propeller  to  turn  one  revolution.  For  simplicity,  the  characteristic  length  L0 
was  chosen  to  be  1  ft  so  the  full  scale  tugboat  and  propeller  grids  (in  ft)  can  be  used  directly  without 
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rescaling.  This  gives  a  Reynolds  number  of  2.647  x  105  based  on  the  characteristic  length  L0  =  1  ft 
when  the  propeller  is  rotating  at  200  rpm.  The  corresponding  Reynolds  number  based  on  the 
propeller  diameter  is  1.488  x  107  based  on  the  propeller  diameter.  The  RANS  equations  are  solved 
with  the  near-wall  two-layer  k-s  turbulence  model. 

The  movie_x*  .  dat  andmovie_q*  .  dat  output  files  were  post-processed  using  the  FieldView 
flow  visualization  software.  Typical  results  include  the  velocity  contours  and  velocity  vector  plots  at 
selected  coordinate  surfaces,  and  pressure  contours  on  the  propeller  blade  and  shroud  surfaces.  The 
velocity  and  pressure  fields  induced  by  the  twin  propellers  were  shown  earlier  in  Figures  19  through 
23  for  this  case.  These  velocity  contours  and  velocity-vector  plots  clearly  illustrate  that  the  right 
propeller  wake  is  strongly  affected  by  the  parallel  pier  wall.  Furthermore,  there  is  a  strong  interaction 
between  the  left  and  right  propellers  with  the  two  ducted  propellers  rotating  in  the  same  rotation.  For 
the  co-rotating  propellers  considered  here,  there  is  a  partial  suppression  of  the  swirling  flow 
momentums  in  the  overlap  region  between  two  propeller  wakes.  This  resulted  in  a  deflection  of  the 
weaker  left  propeller  wake  (away  from  the  pier  wall)  toward  the  sea  bottom,  as  shown  in  Figure  23 
provided  in  the  SSC  Pacific  technical  report. 
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APPENDIX  A 

STRUCTURE  OF  THE  BOUNDARY  CONDITION  INPUT 

The  structure  of  the  boundary  conditions  input  is  best  explained  by  showing  the  pseudo-code  used 
in  the  program  to  read  the  data: 

do  nbk=l , nblocks 

read (LB,  '  (a40)  ')  blocknameO 

!  dummy  read,  blockname  specified  in  inputblk.dat 
read (LB,*)  mb, nf abcs (mb) , ni_dum, n j  dum, nk  dum, mlamp (mb) 

read (LB, * ) 

read (LB,*)  nregions (mb) 

nregs=nregions (mb) 
do  nr=l,nregs 

read (LB, * ) nreg, iedy (mb, nr) , idist (mb, nr) ,  & 
read (LB, * ) mstl (mb, nr) , mndl (mb,  nr)  ,  & 
mst2 (mb, nr ) , mnd2 (mb, nr ) ,  & 
mst3 (mb, nr) , mnd3 (mb, nr) 

end  do 
read (LB, * ) 

read (LB,*)  nLSregions (mb) 

nLSregs=nLSregions (mb) 
do  nLSr=l , nLSregs 

read (LB, * ) nLSreg, nLS (mb, nLSr) , 

mstl  LS (mb, nLSr ) , mndl_LS (mb, nLSr )  ,  & 
mst2_LS (mb, nLSr) , mnd2_LS (mb, nLSr) ,  & 
mst3_LS (mb, nLSr) , mnd3_LS (mb, nLSr) 

end  do 

read (LB, * )  mb,  nfabcs (mb) , nregions  (mb) 
nfbcs=nfabcs (mb) 
do  nf=l,nfbcs 

read (LB, * ) nf a (mb, nf) ,nsec(mb,nf) 
nsect=nsec (mb, nf ) 
do  ns=l,nsect 

read (LB, * ) nu (mb, nf , ns ) , nv (mb, nf , ns ) , nw (mb, nf , ns ) , npr (mb, nf , ns ) ,  & 
ntke (mb, nf , ns) , ntds (mb, nf , ns) , nphiLS (mb, nf , ns) 
read (LB, *) nstl (mb, nf , ns) , nndl (mb, nf ,  ns)  ,  & 
nst2 (mb, nf , ns) , nnd2 (mb, nf , ns) 
end  do 
end  do 

end  do 

read (LB, * )  nbk_prd,  i_prd, j_prd, k_prd 
close (LB) 

Each  process  expects  a  boundary  condition  input  file  and  executes  the  above  given  pseudo-code. 
The  first  read  statement  is  a  dummy  read  of  the  blockname  to  provide  clarity,  while  the 
blockname  specified  earlier  in  the  inputblk .  dat  file  will  be  used  for  block  identification. 
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The  second  read  statement  reads-in  the  local  block  number,  the  number  of  faces  with  boundary 
conditions  for  that  block,  and  the  (i,  j,  k)  dimensions  of  the  block. 

The  third  read  statement  reads-in  the  number  of  regions  needed  to  identify  near-wall  and  outer 
regions  for  the  two-layer  k-s  model.  For  each  region,  we  read  the  iedy  flag,  indicating  whether  the 
region  is  a  near-wall  region  (a  value  of  1)  or  an  outer  region  (a  value  of  2).  If  the  region  is  a  near-wall 
region  then  the  value  of  idi  st  is  the  block  face  number  on  which  the  wall  lies.  We  then  read  the  (i, 
j,  k)  size  of  the  region. 

The  fourth  read  statement  reads-in  the  number  of  regions  needed  for  various  treatments  of  level-set 
function  for  free  surface  flows.  For  each  region,  we  read  the  nLS  flag,  indicating  whether  the  free 
surface  should  be  solved  directly  using  the  advection  equation  for  level-set  function  (a  value  of  1), 
updated  using  zero-gradient  condition  (a  value  of  2,  typically  used  for  the  near-wall  region),  or 
skipped  (a  value  of  3,  for  single-phase  regions  without  air- water  interface).  The  (i,  j,  k)  size  was  then 
specified  for  each  region.  For  the  propeller  wash  study  considered  here,  the  free  surface  wave  effects 
were  negligible  and  the  initial  level-set  function  for  calm  free  surface  was  used  throughout  the  entire 
simulation  by  specifying  nLS  =3. 

After  defining  various  regions  for  the  two-layer  k-s  turbulence  model  and  the  level-set  function, 
we  then  specify  boundary  conditions  for  all  boundary  surfaces  in  the  following  order:  i  =  imin  (Face 
#1),  i  =  imax  (Face  #2),  j  =  jmin  (Face  #3),  j  =jmax  (Face  #4),  k  =  kmin  (Face  #5),  and  k  =  kmax 
(Face  #6).  For  each  face,  we  first  read  the  face  number  and  the  number  of  sections  in  the  face.  For 
each  section  on  a  given  face,  we  then  read  7  boundary  conditions  associated  with  each  of  the  7  field 
variables:  (u,  v,  w,  p,  k,  s,  phiLS),  and  the  surface  limits  on  that  face.  Once  this  is  done  for  all  faces, 
we  proceed  to  read  the  block,  region,  and  surface  data  for  the  next  local  block  on  the  same  processor 
(if  more  than  one  blocks  are  assigned  to  the  same  CPU).  Finally,  we  must  specify  where  the  global 
pressure  datum  is  located.  The  information  is  stored  in  the  variables,  nbk_prd,  i_prd,  j_prd, 
and  k  prd. 

A  list  of  available  boundary  conditions  for  the  velocity  components  and  turbulence  field  variables 
is  given  below: 

•  #1 :  Dirichlet  boundary  condition,  which  is  set  by  the  initial  guess  or  in  the  initial  input 

•  #2:  linear-extrapolation  boundary  condition 

•  #3:  homogeneous  Neumann  or  zero  gradient  boundary  condition 

•  #4:  interior  boundary  condition  for  overset  grids,  interpolation  using  donor  data 

•  #5:  prescribed  boundary  condition,  which  is  specified  by  initial  input  or  updated  in  flow  solver 

•  #6:  moving  surface  boundary  condition,  assign  grid  velocities 

•  #7,  #8:  free  (not  used  at  the  moment) 

•  #9:  moving  surface  boundary  condition,  assign  grid  velocities 

•  #10:  free  (not  used  at  the  moment) 

•  #11:  branch  cut  in  lower  index,  average  across  branch  cut 

•  #12:  branch  cut  in  higher  index,  average  across  branch  cut 

•  #13:  collapse-to-axis  in  lower  index,  average  in  circumferential  direction 

•  #14:  collapse-to-axis  in  higher  index,  average  in  circumferential  direction 


The  following  is  a  list  of  available  boundary  conditions  for  the  pressure: 

•  #1 :  free  (not  used  at  the  moment) 

•  #2:  linear-extrapolation  boundary  condition 


A-2 


•  #3:  homogeneous  Neumann  or  zero  gradient  boundary  condition 

•  #4:  interior  boundary  condition  for  overset  grids,  interpolation  using  donor  data 

•  #5:  prescribed  boundary  condition,  which  is  specified  by  initial  input  or  updated  in  flow  solver 

•  #6:  free  (not  used  at  the  moment) 

•  #7:  compute  pressure  consistently,  using  conservation  of  mass 

•  #8,  #9,  #10:  free  (not  used  at  the  moment) 

•  #11:  branch  cut  in  lower  index,  average  across  branch  cut 

•  #12:  branch  cut  in  higher  index,  average  across  branch  cut 

•  #13:  collapse-to-axis  in  lower  index,  average  in  circumferential  direction 

•  #14:  collapse-to-axis  in  higher  index,  average  in  circumferential  direction 

Typical  boundary  conditions  for  a  stationary  wall  are  either  of  the  following: 


1112113  or 
1113  113 


In  the  first  one,  pressure  is  linearly  extrapolated  at  the  wall,  and  in  the  second  it  is  computed 
consistently  at  the  wall  using  conservation  of  mass  at  the  boundary  itself.  For  highly  skewed  meshes 
in  the  near  wall  region,  linear  extrapolation  is  more  stable. 

At  a  free-stream  inflow  the  following  are  valid  options: 

5  5  5  2  1  1  or 

5  5  5  3  1  1 


In  the  first  one,  a  zero  pressure  gradient  is  enforced,  and  in  the  second  pressure  is  computed 
consistently  at  the  boundary  using  conservation  of  mass.  The  second  option,  where  pressure  is 
computed  consistently,  is  also  valid  at  the  inflow  of  a  channel  -  where  a  pressure  drop  is  present.  The 
user  can  appreciate  the  versatility  of  the  consistent  pressure  boundary  condition,  as  it  applies  to 
virtually  any  situation  where  velocities  are  prescribed. 

Similarly,  at  an  outflow,  the  following  are  valid  options: 

2  2  2  2  2  2  or 

2  2  2  7  2  2 


For  problems  involving  free  surface,  the  available  boundary  conditions  for  the  level-set  function 
are  listed  in  the  following: 

•  #1 :  free  (not  used  at  the  moment) 

•  #2 :  linear-extrapolation  boundary  condition 

•  #3:  homogeneous  Neumann  or  zero  gradient  boundary  condition 

•  #4:  interior  boundary  condition  for  overset  grids,  interpolation  using  donor  data 

•  #5:  prescribed  boundary  condition,  which  is  specified  by  initial  input  or  updated  in  flow  solver 

•  #6,  #7,  #8,  #9,  #10:  free  (not  used  at  the  moment) 

•  #11:  branch  cut  in  lower  index,  average  across  branch  cut 

•  #12:  branch  cut  in  higher  index,  average  across  branch  cut 
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#13:  collapse-to-axis  in  lower  index,  average  in  circumferential  direction 
#14:  collapse-to-axis  in  higher  index,  average  in  circumferential  direction 


APPENDIX  B 

COSMIC  INPUT  FILE  FOR  DDG-51  AND  P4876 
PROPELLER  WASH  STUDY 


!  example  input  for  DDG-51  Ship  and  P4876  propeller  wash  study 

!  global  parameters 

<global> 

fringe  =  1, 
quality  =  0.01, 
nquality  =  4, 
eps  =  0.001, 

</ global> 


!  grid  block  definition 

<block  name  =  "propellerO 1 "> 

linking  grid  list  =  <"propeller01" , "shaftOl", "shaft02", "basinOl", 
"basin02", "propeller05" , "propeller02 " , "ocean03"> 

</block> 

<block  name  =  "propeller02 "> 

linking  grid  list  =  <"propeller02 " , "shaf tO 1 " , "shaf t02 " , "basinO 1 " , 
"basin02", "propellerO 1 " , "propeller03" , "ocean03"> 

</block> 

<block  name  =  "propeller03"> 

linking  grid_list  =  <"propeller03" , "shaf tOl ", "shaf t02 ", "basinOl "  , 
"basin02", "propeller02 " , "propeller04", "ocean03"> 

</block> 

<block  name  =  "propeller04"> 

linking  grid  list  =  <"propeller04" , "shaftOl", "shaft02", "basinOl", 
"basin02", "propeller03" , "propeller05" , "ocean03"> 

</block> 

<block  name  =  "propeller05"> 

linking  grid  list  =  <"propeller05" , "shaftOl", "shaft02",  "basinOl", 
"basin02", "propeller04", "propellerO 1 " , "ocean03"> 

</block> 

<block  name  =  "shaft01"> 

linking  grid  list  =  <"shaf tO 1 ", "shaf t02 ", "shaf t03"  ,  "basinO 1 "  , 
"basin02", "propellerO 1 " , "propeller02 " , 

"propeller03" , "propeller04", "propeller05" , 

"ocean02", "ocean03"> 

</block> 

<block  name  =  "shaft02"> 

linking  grid  list  =  <"shaf t02 ", "shaf tO 1 ", "shaf t03" , "basinO 1 "  , 
"basin02", "propellerO 1 " , "propeller02 " , 

"propeller03" , "propeller04", "propeller05" , 

"ocean02", "ocean03"> 

</block> 

<block  name  =  "shaft03"> 

linking  grid_list  =  <"shaft03" , "shaf tOl ", "shaf t02 ", "basinOl "  , 
"shipOl", "ocean02", "ocean03"> 

</block> 

<block  name  =  "ship01"> 

linking  grid_list  =  <"basin0 1 " , "basin02 " , "basin03"  ,  "oceanO 1 "  , 


B-l 


"ocean02 " , "ocean03" , " shaftOl " , " shaft 02 " , 

"shaf t03"> 

</block> 

<block  name  =  "basin01"> 

linking  grid_list  =  <"basin02 ", "shaf tO 1 ", "shaf t02 ", "shaf t03" , 
"oceanOl", "ocean02", "ocean03", 

"propellerO 1 " , "propeller02 " , "propeller03" , 
"propeller04" , "propeller05" , "ship01"> 

</block> 

<block  name  =  "basin02"> 

linking  grid_list  =  <"basin0 1 ", "basin03" , "shaf tO 1 ", "shaf t02 " , 
"shaf tO 3" , "oceanO 1 " , "ocean02 " , "ocean03" , 
"propellerO 1 " , "propeller02 " , "propeller03" , 
"propeller04" , "propeller05" , "ship01"> 

</block> 

<block  name  =  "basin03"> 

linking  grid  list  =  <"basin02 " , "oceanOl " , "ocean03"> 

</block> 

<block  name  =  "ocean01"> 

linking  grid  list  =  <"basin01 "  ,  "basin02 " , "basin03" , 

"ocean02", "ocean03"> 

</block> 

<block  name  =  "ocean02"> 

linking  grid  list  =  <"ocean0 1 " , "ocean03" , "shipOl " , "basinOl "  , 
"shaftOl", "shaft02", "shaft03"> 

</block> 

<block  name  =  "ocean03"> 

linking  grid  list  =  <"ocean0 1 " , "ocean02 " , "basinO 1 " , "basin02 " , 
"basin03", "shipOl", "shaft03", "shaft02", 

" shaftOl "> 

</block> 

<block  name  =  "phantom01"> 

linking  grid  list  =  <"phantom01"> 

</block> 

<block  name  =  "phantom02"> 

linking  grid_list  =  <"phantom02 "> 

</block> 

<block  name  =  "phantom03"> 

linking  grid  list  =  <"phantom03"> 

</block> 

<block  name  =  "phantom04"> 

linking  grid  list  =  <"phantom04"> 

</block> 

<block  name  =  "phantom05"> 

linking  grid_list  =  <"phantom05"> 

</block> 

<block  name  =  "phantom06"> 

linking  grid  list  =  <"phantom06"> 

</block> 

<block  name  =  "phantom07"> 

linking  grid_list  =  <"phantom07"> 

</block> 


!  hole  boundary  definition 

<boundary  name  =  "phantomOl  hole  boundary"> 
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parent  grid  =  "phantomOl", 

hole  cutting  list  =  <"shaf tOl " , "shaf t02 " , "basinOl " , "basin02 "> 
</boundary> 

<boundary  name  =  "phantom02  hole  boundary"> 
parent  grid  =  "phantom02", 

hole  cutting_list  =  <"shaft01", "shaft02", "basinOl", "basin02"> 
</boundary> 

<boundary  name  =  "phantom03  hole  boundary"> 
parent  grid  =  "phantom03", 

hole  cutting_list  =  <"shaft01",  "shaft02", "basinOl", "basin02"> 
</boundary> 

<boundary  name  =  "phantom04  hole  boundary"> 
parent  grid  =  "phantom04", 

hole  cutting  list  =  <"shaf tO 1 " ,  "shaf t02 " , "basinO 1 " , "basin02 "> 
</boundary> 

<boundary  name  =  "phantom05  hole  boundary"> 
parent  grid  =  "phantom05", 

hole  cutting_list  =  <"shaft01", "shaft02", "basinOl", "basin02"> 
</boundary> 

<boundary  name  =  "phantom06  hole  boundary"> 
parent  grid  =  "phantom06", 
hole_cutting  list  =  <"ship01"> 

</boundary> 

<boundary  name  =  "phantom07  hole  boundary"> 
parent  grid  =  "phantom07", 
hole^cutting  list  =  <"ship01", "ocean03"> 

</boundary> 

<boundary  name  =  "shaftOl  hole  boundary"> 
parent  grid  =  "shaftOl" 

hole_cutting  list  =  <"basin01", "basin02"> 

</boundary> 

<boundary  name  =  "shaft03  hole  boundary"> 
parent^grid  =  "shaft03" 
hole_cutting  list  =  <"basin01"> 

</boundary> 

<boundary  name  =  "shipOl  hole  boundary"> 
parent  grid  =  "shipOl", 

hole_cutting  list  =  <"shaft03"  ,  "basinOl ", "basin02 " , 
"ocean02", "ocean03"> 

</boundary> 


!  hole  surface  definitions 

<surface  name  =  "phantomOl  hole  boundary"> 
ijk^range  =  1,  1,  1,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk". 
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</surf ace> 

<surface  name  =  "phantomOl  hole  boundary"> 
ijk_range  =  3,  3,  1,  41,  1,  61, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantomOl  hole  boundary"> 
ijk_range  =  1,  3,  1,  1,  1,  61, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantomOl  hole  boundary"> 
ijk_range  =  1,  3,  41,  41,  1,  61, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantomOl  hole  boundary"> 
ijk^range  =  1,  3,  1,  41,  61,  61, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom02  hole  boundary"> 
ijk_range  =  1,  1,  1,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom02  hole  boundary"> 
ijk_range  =  3,  3,  1,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom02  hole  boundary"> 
ijk^range  =  1,  3,  1,  1,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom02  hole  boundary"> 
ijk_range  =  1,  3,  41,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom02  hole  boundary"> 
ijk^range  =  1,  3,  1,  41,  61,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom03  hole  boundary"> 
ijk^range  =  1,  1,  1,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom03  hole  boundary"> 
ijk^range  =  3,  3,  1,  41,  1,  61, 
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boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom03  hole  boundary"> 
ijk^range  =  1,  3,  1,  1,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom03  hole  boundary"> 
ijk_range  =  1,  3,  41,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom03  hole  boundary"> 
ijk_range  =  1,  3,  1,  41,  61,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom04  hole  boundary"> 
ijk_range  =  1,  1,  1,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom04  hole  boundary"> 
ijk_range  =  3,  3,  1,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom04  hole  boundary"> 
ijk_range  =  1,  3,  1,  1,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom04  hole  boundary"> 
ijk^range  =  1,  3,  41,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom04  hole  boundary"> 
ijk_range  =  1,  3,  1,  41,  61,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom05  hole  boundary"> 
ijk^range  =  1,  1,  1,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom05  hole  boundary"> 
ijk^range  =  3,  3,  1,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 
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<surface  name  =  "phantom05  hole  boundary"> 
ijk^range  =  1,  3,  1,  1,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom05  hole  boundary"> 
ijk^range  =  1,  3,  41,  41,  1,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom05  hole  boundary"> 
ijk^range  =  1,  3,  1,  41,  61,  61, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom06  hole  boundary"> 
ijk_range  =  1,  1,  1,  2,  1,  2, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom06  hole  boundary"> 
ijk_range  =  2,  2,  1,  2,  1,  2, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom06  hole  boundary"> 
ijk_range  =  1,  2,  1,  1,  1,  2, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom06  hole  boundary"> 
ijk_range  =  1,  2,  2,  2,  1,  2, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom06  hole  boundary"> 
ijk_range  =  1,  2,  1,  2,  1,  1, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom06  hole  boundary"> 
ijk_range  =  1,  2,  1,  2,  2,  2, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom07  hole  boundary"> 
ijk^range  =  1,  1,  1,  2,  1,  2, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom07  hole  boundary"> 
ijk_range  =  2,  2,  1,  2,  1,  2, 

boundary  condition  ="cut". 
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surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom07  hole  boundary"> 
ijk_range  =  1,  2,  1,  1,  1,  2, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom07  hole  boundary"> 
ijk_range  =  1,  2,  2,  2,  1,  2, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom07  hole  boundary"> 
ijk_range  =  1,  2,  1,  2,  1,  1, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom07  hole  boundary"> 
ijk_range  =  1,  2,  1,  2,  2,  2, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "shaftOl  hole  boundary"> 

ijk_range  =  1,  38,  3,  3,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "shaft03  hole  boundary"> 

ijk_range  =  6,  65,  3,  3,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "shipOl  hole  boundary"> 

ijk  range  =  1,  121,  21,  21,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


!  outer  boundary  definition 

<boundary  name  =  "propellerOl  outer  boundary"> 
parent  grid  =  "propellerOl", 

</boundary> 

<boundary  name  =  "propeller02  outer  boundary"> 
parent  grid  =  "propeller02 " , 

</boundary> 

<boundary  name  =  "propeller03  outer  boundary"> 
parent  grid  =  "propeller03" , 

</boundary> 
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<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent_grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 
<boundary  name  =  ' 
parent  grid 
</boundary> 


propeller04  outer  boundary"> 
=  "propeller04" , 

propeller05  outer  boundary"> 
=  "propeller05" , 

shaftOl  outer  boundary"> 

=  "shaftOl" , 

shaft02  outer  boundary"> 

=  "shaft02", 

shaft03  outer  boundary"> 

=  "shaft03". 


shipOl  outer  boundary"> 

=  "shipOl", 

basinOl  outer  boundary"> 
=  "basinOl", 

basin02  outer  boundary"> 
=  "basin02", 

basin03  outer  boundary"> 
=  "basin03". 


oceanOl  outer  boundary"> 
=  "oceanOl", 


ocean02  outer  boundary"> 
=  "ocean02". 


ocean03  outer  boundary"> 
=  "ocean03". 


!  outer  boundary  surface  definition 

<surface  name  =  "propellerOl  outer  boundary"> 
ijk  range  =  1,  1,  1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propellerOl" 
donor_ijk  range  =  61,  61,  1,  41,  1,  41, 

</surf ace> 

<surface  name  =  "propellerOl  outer  boundary"> 
ijk_range  =  62,  62,1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propellerOl" 
donor_ijk  range  =  2,2,1,  41,1,  41, 

</surf ace> 

<surface  name  =  "propellerOl  outer  boundary"> 
ijk_range  =  1,  62,  41,  41,  2,  40, 

</surf ace> 
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<surface  name  =  "propeller02  outer  boundary"> 
ijk  range  =  1,  1,  1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller02 " 
donor_ijk  range  =  61,  61,  1,  41,  1,  41, 

</surf ace> 

<surface  name  =  "propeller02  outer  boundary"> 
ijk_range  =62,  62,1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller02 " 
donor_ijk  range  =  2,2,1,  41,1,  41, 

</surf ace> 

<surface  name  =  "propeller02  outer  boundary"> 
ijk_range  =  1,  62,  41,  41,  2,  40, 

</surf ace> 


<surface  name  =  "propeller03  outer  boundary"> 
ijk  range  =  1,  1,  1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller03" 
donor_ijk  range  =  61,  61,  1,  41,  1,  41, 

</surf ace> 

<surface  name  =  "propeller03  outer  boundary"> 
ijk_range  =62,  62,1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller03" 
donor_ijk  range  =  2,2,1,  41,1,  41, 

</surf ace> 

<surface  name  =  "propeller03  outer  boundary"> 
ijk_range  =  1,  62,  41,  41,  2,  40, 

</surf ace> 


<surface  name  =  "propeller04  outer  boundary"> 
ijk  range  =  1,  1,  1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller04" 
donor_ijk  range  =  61,  61,  1,  41,  1,  41, 

</surf ace> 

<surface  name  =  "propeller04  outer  boundary"> 
ijk_range  =62,  62,1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller04" 
donor_ijk  range  =  2,2,1,  41,1,  41, 

</surf ace> 

<surface  name  =  "propeller04  outer  boundary"> 
ijk_range  =  1,  62,  41,  41,  2,  40, 

</surf ace> 


<surface  name  =  "propeller05  outer  boundary"> 
ijk  range  =  1,  1,  1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller05" 
donor_ijk  range  =  61,  61,  1,  41,  1,  41, 

</surf ace> 
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<surface  name  =  "propeller05  outer  boundary"> 
ijk_range  =  62,  62,1,  41,  1,  41, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller05" 
donor_ijk  range  =  2,2,1,  41,1,  41, 

</surf ace> 

<surface  name  =  "propeller05  outer  boundary"> 
ijk_range  =  1,  62,  41,  41,  2,  40, 

</surf ace> 


<surface  name  =  "shaftOl  outer  boundary"> 
ijk  range  =  1,  1,  1,  21,  1,  122, 

</surf ace> 

<surface  name  =  "shaftOl  outer  boundary"> 
ijk  range  =  1,  38,  21,  21,  1,  122, 

</surf ace> 

<surface  name  =  "shaftOl  outer  boundary"> 
ijk_range  =  1,  38,  1,  21,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "shaftOl" 

donor_ijk  range  =  1,  38,  1,  21,  121,  121, 

</surf ace> 

<surface  name  =  "shaftOl  outer  boundary"> 
ijk  range  =  1,  38,  1,  21,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "shaftOl" 

donor_ijk  range  =  1,  38,  1,  21,  2,  2, 

</surf ace> 


<surface  name  =  "shaft02  outer  boundary"> 
ijk  range  =  1,  1,  1,  31,  1,  122, 

</surf ace> 

<surface  name  =  "shaft02  outer  boundary"> 
ijk_range  =28,  28,  1,  31,  1,  122, 

</surf ace> 

<surface  name  =  "shaft02  outer  boundary"> 
ijk  range  =  1,  28,  1,  1,  1,  122, 

</surf ace> 

<surface  name  =  "shaft02  outer  boundary"> 
ijk_range  =  1,  28,  31,  31,  1,  122, 

</surf ace> 

<surface  name  =  "shaft02  outer  boundary"> 
ijk^range  =  1,  28,  1,  31,  1,  1, 

boundary^condition  =  "periodic", 
donor^grid  =  "shaft02" 

donor^ijk  range  =  1,  28,  1,  31,  121,  121, 

</surf ace> 

<surface  name  =  "shaft02  outer  boundary"> 
ijk_range  =  1,  28,  1,  31,  122,  122, 

boundary^condition  =  "periodic", 
donor^grid  =  "shaft02" 

donor^ijk  range  =  1,  28,  1,  31,  2,  2, 

</surf ace> 
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<surface  name  =  "shaft03  outer  boundary"> 
ijk_range  =65,  65,  1,  21,  1,  122, 

</surf ace> 

<surface  name  =  "shaft03  outer  boundary"> 
ijk  range  =  2,  65,  21,  21,  1,  122, 

</surf ace> 

<surface  name  =  "shaft03  outer  boundary"> 
ijk_range  =  1,  65,  1,  21,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "shaft03" 

donor_ijk  range  =  1,  65,  1,  21,  121,  121, 

</surf ace> 

<surface  name  =  "shaft03  outer  boundary"> 
ijk  range  =  1,  65,  1,  21,  122,  122, 
boundary  condition  =  "periodic", 
donor_grid  =  "shaft03" 

donor_ijk  range  =  1,  65,  1,  21,  2,  2, 

</surf ace> 

<surface  name  =  "shaft03  outer  boundary"> 
ijk_range  =  1,  65,  1,  1,  1,  122, 

boundary  condition  =  "body", 

</surf ace> 


<surface  name  =  "shipOl  outer  boundary"> 

ijk_range  =  1,  121,  35,  35,  1,  41, 

</surf ace> 


<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 


"basinOl  outer  boundary"> 
=  1,  1,  1,  81,  1,  77, 

"basinOl  outer  boundary"> 
=  34,  34,  1,  81,  1,  77, 

"basinOl  outer  boundary"> 
=  1,  34,  81,  81,  1,  77, 

"basinOl  outer  boundary"> 
=  1,  34,  1,  81,  1,  1, 

"basinOl  outer  boundary"> 
=  1,  34,  1,  81,  77,  77, 


<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 


"basin02  outer  boundary"> 
=  1,  1,  1,  81,  1,  77, 

"basin02  outer  boundary"> 
=  34,  34,  1,  81,  1,  77, 

"basin02  outer  boundary"> 
=  1,  34,  81,  81,  1,  77, 

"basin02  outer  boundary"> 
=  1,  34,  1,  81,  1,  1, 
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</surf ace> 

<surface  name  =  "basin02  outer  boundary"> 
ijk_range  =  1,  34,  1,  81,  77,  77, 

</surf ace> 


<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 


"basin03  outer  boundary"> 
=  1,  1,  1,  81,  1,  77, 

"basin03  outer  boundary"> 
=  34,  34,  1,  81,  1,  77, 

"basin03  outer  boundary"> 
=  1,  34,  81,  81,  1,  77, 

"basin03  outer  boundary"> 
=  1,  34,  1,  81,  1,  1, 

"basin03  outer  boundary"> 
=  1,  34,  1,  81,  77,  77, 


<surface  name  =  "oceanOl  outer  boundary"> 
ijk_range  =  1,  152,  1,  65,  21,  21, 

</surf ace> 


<surface  name  =  "ocean02  outer  boundary"> 
ijk_range  =77,  77,  1,  65,  1,  42, 

</surf ace> 

<surface  name  =  "ocean02  outer  boundary"> 
ijk_range  =  1,  77,  1,  65,  1,  1, 

</surf ace> 


<surface  name  =  "ocean03  outer  boundary"> 
ijk_range  =  1,  1,  1,  65,  1,  42, 

</surf ace> 

<surface  name  =  "ocean03  outer  boundary"> 
ijk_range  =  1,  77,  1,  65,  1,  1, 

</surf ace> 
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APPENDIX  C 

COSMIC  INPUT  FILE  FOR  TUGBOAT 
AND  DUCTED-PROPELLER  WASH  STUDY 


!  example  input  for  tugboat  and  ducted  propellers 

!  global  parameters 

<global> 

fringe  =  1, 
quality  =  0.01, 
nquality  =  4, 
eps  =  0.001, 

</ global> 


!  grid  block  definition 
<block  name  =  "propellerO la"> 

linking  grid  list  =  <"propeller0 la" , "tipOla" , "propeller04a" , 
"propeller02a" , "duct02a", "duct03a"> 

</block> 

<block  name  =  "propeller02a"> 

linking  grid  list  =  <"propeller02a" , "tip02a" , "propellerO la" , 
"propeller03a" , "duct02a", "duct03a"> 

</block> 

<block  name  =  "propeller03a"> 

linking  grid_list  =  <"propeller03a" , "tip03a" , "propeller02a"  , 
"propeller04a" , "duct02a", "duct03a"> 

</block> 

<block  name  =  "propeller04a"> 

linking  grid  list  =  <"propeller04a" , "tip04a" , "propeller03a"  , 
"propellerO la" , "duct02a", "duct03a"> 

</block> 

<block  name  =  "tip01a"> 

linking  grid_list  =  <"propeller0 la" , "duct03a"> 

</block> 

<block  name  =  "tip02a"> 

linking  grid  list  =  <"propeller02a" , "duct03a"> 

</block> 

<block  name  =  "tip03a"> 

linking  grid_list  =  <"propeller03a" , "duct03a"> 

</block> 

<block  name  =  "tip04a"> 

linking  grid  list  =  <"propeller04a" , "duct03a"> 

</block> 

<block  name  =  "duct01a"> 

linking  grid_list  =  <"duct01a" , "duct02a" , "duct03a"  ,  "duct04a"  , 
"bargeOl", "ocean03", "ocean04"> 

</block> 

<block  name  =  "duct02a"> 

linking  grid  list  =  <"duct02a" , "ductOla" , "duct03a"  ,  "duct05a"  , 
"wakeOla", "propellerO la" , "propeller02a"  , 
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"propeller03a" , "propeller04a"> 

</block> 

<block  name  =  "duct03a"> 

linking  grid_list  =  <"duct03a" , "ductOla" , "duct02a" , "duct05a" , 
"propellerO la" , "propeller02a" , 

"propeller03a" , "propeller04a" , "tipOla", 

"tip02a", "tip03a", "tip04a", "barge01"> 

</block> 

<block  name  =  "duct04a"> 

linking  grid  list  =  <"duct04a" , "ductOla" , "duct05a"  ,  "wakeOla"  , 
"bargeOl", "ocean04"> 

</block> 

<block  name  =  "duct05a"> 

linking  grid  list  =  <"duct05a" , "duct02a" , "duct03a"  ,  "duct04a"  , 
"wakeOla", "bargeOl", "ocean04"> 

</block> 

<block  name  =  "wake01a"> 

linking  grid  list  =  <"wake01a" , "duct02a" , "duct03a"  ,  "duct04a"  , 
"duct05a", "bargeOl", "ocean04", "ocean05"> 

</block> 

<block  name  =  "propellerO lb"> 

linking  grid  list  =  <"propellerO lb" , "tipO lb" , "propeller04b" , 
"propeller02b" , "duct02b", "duct03b"> 

</block> 

<block  name  =  "propeller02b"> 

linking  grid_list  =  <"propeller02b" , "tip02b" , "propellerO lb" , 
"propeller03b" , "duct02b", "duct03b"> 

</block> 

<block  name  =  "propeller03b"> 

linking  grid  list  =  <"propeller03b" , "tip03b" , "propeller02b" , 
"propeller04b" , "duct02b", "duct03b"> 

</block> 

<block  name  =  "propeller04b"> 

linking  grid  list  =  <"propeller04b" , "tip04b" , "propeller03b" , 
"propellerO lb" , "duct02b", "duct03b"> 

</block> 

<block  name  =  "tip01b"> 

linking  grid_list  =  <"propellerO lb" , "duct03b"> 

</block> 

<block  name  =  "tip02b"> 

linking  grid  list  =  <"propeller02b" , "duct03b"> 

</block> 

<block  name  =  "tip03b"> 

linking  grid_list  =  <"propeller03b" , "duct03b"> 

</block> 

<block  name  =  "tip04b"> 

linking  grid_list  =  <"propeller04b" , "duct03b"> 

</block> 

<block  name  =  "duct01b"> 

linking  grid_list  =  <"duct01b" , "duct02b" , "duct03b" , "duct04b" , 
"bargeOl", "ocean03", "ocean04"> 

</block> 

<block  name  =  "duct02b"> 

linking  grid  list  =  <"duct02b" , "ductOlb" , "duct03b"  ,  "duct05b"  , 
"wakeOlb", "propellerO lb" , "propeller02b" , 

"propeller 03b" , "propeller04b"> 

</block> 
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<block  name  =  "duct03b"> 

linking  grid  list  =  <"duct03b" , "ductOlb" , "duct02b" , "duct05b" , 
"propellerO lb" , "propeller 02b" , 

"propeller03b" , "propeller04b" , "tipOlb", 

"tip02b", "tip03b", "tip04b", "barge01"> 

</block> 

<block  name  =  "duct04b"> 

linking  grid  list  =  <"duct04b" , "ductOlb" , "duct05b" , "wakeOlb" , 
"bargeOl", "ocean04"> 

</block> 

<block  name  =  "duct05b"> 

linking  grid_list  =  <"duct05b" , "duct02b" , "duct03b" , "duct04b" , 
"wakeOlb", "bargeOl", "ocean04"> 

</block> 

<block  name  =  "wake01b"> 

linking  grid  list  =  <"wake01b" , "duct02b" , "duct03b" , "duct04b"  , 
"duct05b" , "bargeOl" , "ocean04" , "ocean05"> 

</block> 

<block  name  =  "ocean01"> 

linking  grid^list  =  <"ocean02 " , "oceanl7 " , "oceanl8 "  ,  "bargeO 1 "> 

</block> 

<block  name  =  "ocean02"> 

linking  grid_list  =  <"oceanO 1 " , "ocean03" , "oceanl7 "  ,  "oceanl 8 "  , 
"bargeOl "> 

</block> 

<block  name  =  "ocean03"> 

linking  grid  list  =  <"ocean02 ", "ocean04 ", "oceanl7 ", "oceanl 8 "  , 
"ductOla", "ductOlb", "duct04a", "duct04b", 

"bargeOl "> 

</block> 

<block  name  =  "ocean04"> 

linking  grid  list  =  <"ocean03" , "ocean05" , "oceanl7 ", "oceanl 8 " , 
"wakeOla", "wakeOlb", "ductOla", "ductOlb", 

"duct04a", "duct04b", "duct05a", "duct05b", 

"bargeOl "> 

</block> 

<block  name  =  "ocean05"> 

linking  grid  list  =  <"ocean04 " , "oceanO 6" , "oceanl7 " , "oceanl 8 "  , 
"wakeOla", "wakeOlb", "barge01"> 

</block> 

<block  name  =  "ocean06"> 

linking  grid  list  =  <"ocean05" , "ocean07 ", "oceanl7 ",  "oceanl 8 "  , 
"wakeOla", "wakeOlb", "barge01"> 

</block> 

<block  name  =  "ocean07"> 

linking  grid_list  =  <"oceanO 6" , "oceanO 8 " , "oceanl7 "  ,  "oceanl 8 "  , 
"bargeOl "> 

</block> 

<block  name  =  "ocean08"> 

linking  grid_list  =  <"ocean07 ", "oceanO 9" , "oceanl7 ",  "oceanl 8 "  , 
"bargeOl "> 

</block> 

<block  name  =  "ocean09"> 

linking  grid_list  =  <"oceanO 8 " , "oceanlO " , "oceanl7 "  ,  "oceanl8 "  , 
"bargeOl "> 

</block> 

<block  name  =  "oceanlO"> 
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oceanl 1 


oceanl7 


oceanl 8 


linking  grid_list  =  <"ocean09"," 

"bargeO 1 "> 

</block> 

<block  name  =  "oceanll"> 

linking  grid  list  =  <"oceanlO ", "oceanl2 ", "oceanl7 ", "oceanl 8 "> 
</block> 

<block  name  =  "oceanl2"> 

linking  grid_list  =  <"oceanl 1 ", "oceanl3" , "oceanl7 ", "oceanl 8 "> 
</block> 

<block  name  =  "oceanl3"> 

linking  grid  list  =  <"oceanl2 ", "oceanl4 ", "oceanl7 ", "oceanl 8 "> 
</block> 

<block  name  =  "oceanl4"> 

linking  grid_list  =  <"oceanl3" , "oceanl5" , "oceanl7 ",  "oceanl 8 "> 
</block> 

<block  name  =  "oceanl5"> 

linking  grid_list  =  <"oceanl4 ", "oceanl 6" , "oceanl7 ", "oceanl 8 "> 
</block> 

<block  name  =  "oceanl6"> 

linking  grid  list  =  <"oceanl5" , "oceanl7 " , "oceanl8 "> 

</block> 

<block  name  =  "oceanl7"> 

linking  grid  list  =  <"oceanO 1 " , "ocean02 " , "ocean03" , "ocean04 " , 
"ocean05", "ocean06", "ocean07", "ocean08", 

"oceanO  9" , " oceanl 0 " , "oceanl 1 " , "ocean 12 " , 

"oceanl3", "oceanl4", "oceanl5", "oceanl6", 

"oceanl8", "barge01"> 

</block> 

<block  name  =  "oceanl8"> 

linking  grid_list  =  <"oceanO 1 " , "ocean02 " , "ocean03" , "ocean04 "  , 
"oceanO 5" , "oceanO  6" , "oceanO 7 " , "oceanO 8 " , 

"oceanO  9" , "ocean 10 " , "oceanl 1 " , "ocean 12 " , 

"oceanl3", "oceanl4", "oceanl5", "oceanl6", 

"oceanl7", "barge01"> 

</block> 

<block  name  =  "barge01"> 

linking  grid  list  =  <"duct01a" , "ductOlb" , "duct04a" , "duct04b" , 
"duct05a", "duct05b", "wakeOla", "wakeOlb", 

"oceanO 1 " , "ocean 02 " , "oceanO 3 " , "ocean 04 " , 

"oceanO 5" , "oceanO  6" , "oceanO 7 " , "oceanO  8 " , 

"oceanO  9" , "ocean 10 " , "ocean 17 " , "oceanl 8 "> 

</block> 

<block  name  =  "phantomO la"> 

linking  grid_list  =  <"phantom0 la"> 

</block> 

<block  name  =  "phantom02a"> 

linking  grid  list  =  <"phantom02a"> 

</block> 

<block  name  =  "phantom03a"> 

linking  grid  list  =  <"phantom03a"> 

</block> 

<block  name  =  "phantom04a"> 

linking  grid  list  =  <"phantom04a"> 

</block> 

<block  name  =  "phantomO lb"> 

linking  grid_list  =  <"phantom01b"> 

</block> 
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<block  name  =  "phantom02b"> 

linking  grid  list  =  <"phantom02b"> 
</block> 

<block  name  =  "phantom03b"> 

linking  grid  list  =  <"phantom03b"> 
</block> 

<block  name  =  "phantom04b"> 

linking  grid  list  =  <"phantom04b"> 
</block> 

<block  name  =  "phantom05"> 

linking  grid_list  =  <"phantom05"> 
</block> 


!  hole  boundary  definition 

<boundary  name  =  "phantomOla  hole  boundary"> 
parent  grid  =  "phantomOla", 
hole  cutting_list  =  <"duct02a" , "duct03a"> 
</boundary> 

<boundary  name  =  "phantom02a  hole  boundary"> 
parent  grid  =  "phantom02a" , 
hole  cutting_list  =  <"duct02a" , "duct03a"> 
</boundary> 

<boundary  name  =  "phantom03a  hole  boundary"> 
parent  grid  =  "phantom03a" , 
hole  cutting  list  =  <"duct02a" , "duct03a"> 
</boundary> 

<boundary  name  =  "phantom04a  hole  boundary"> 
parent  grid  =  "phantom04a" , 
hole  cutting^list  =  <"duct02a" , "duct03a"> 
</boundary> 

<boundary  name  =  "phantomOlb  hole  boundary"> 
parent  grid  =  "phantomOlb", 
hole  cutting  list  =  <"duct02b" , "duct03b"> 
</boundary> 

<boundary  name  =  "phantom02b  hole  boundary"> 
parent  grid  =  "phantom02b" , 
hole  cutting^list  =  <"duct02b" , "duct03b"> 
</boundary> 

<boundary  name  =  "phantom03b  hole  boundary"> 
parent  grid  =  "phantom03b" , 
hole  cutting_list  =  <"duct02b" , "duct03b"> 
</boundary> 

<boundary  name  =  "phantom04b  hole  boundary"> 
parent  grid  =  "phantom04b" , 
hole  cutting^list  =  <"duct02b" , "duct03b"> 
</boundary> 

<boundary  name  =  "phantom05  hole  boundary"> 
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parent  grid  =  "phantom05", 

hole_cutting_list  =  <"oceanO 1 " , "ocean02 " , "ocean03" , 
"ocean04", "ocean05", "ocean06", 

"ocean07", "ocean08", "ocean09", 

"oceanlO"> 

</boundary> 

<boundary  name  =  "ductOla  hole  boundary" > 
parent  grid  =  "ductOla", 

hole  cutting  list  =  <"ocean03" ,  "ocean04 " , "barged "> 
</boundary> 


<boundary  name  =  "duct04a  hole  boundary"> 
parent  grid  =  "duct04a", 

hole  cutting  list  =  <"wake01a" , "ocean04 " , "ocean03" , 
"bargeO 1 "> 

</boundary> 

<boundary  name  =  "duct05a  hole  boundary"> 
parent_grid  =  "duct05a", 

hole  cutting_list  =  <"wake01a"  ,  "ocean04" , "barge01"> 
</boundary> 

<boundary  name  =  "wakeOla  hole  boundary"> 
parent  grid  =  "wakeOla", 

hole  cutting_list  =  <"ocean04",  "ocean05", "barge01"> 
</boundary> 

<boundary  name  =  "ductOlb  hole  boundary" > 
parent  grid  =  "ductOlb", 

hole  cutting_list  =  <"ocean03" , "ocean04 " , "barged "> 
</boundary> 


<boundary  name  =  "duct04b  hole  boundary"> 
parent  grid  =  "duct04b", 

hole  cutting_list  =  <"wakeO lb" , "ocean04 " , "ocean03" , 
"bargeO 1 "> 

</boundary> 

<boundary  name  =  "duct05b  hole  boundary"> 
parent_grid  =  "duct05b", 

hole  cutting_list  =  <"wake01b"  ,  "ocean04" , "barge01"> 
</boundary> 

<boundary  name  =  "wakeOlb  hole  boundary" > 
parent  grid  =  "wakeOlb", 

hole  cutting^list  =  <"ocean04 ", "ocean05" , "bargeO 1 "> 
</boundary> 

<boundary  name  =  "barged  hole  boundary"> 
parent  grid  =  "bargeOl", 

hole_cutting_list  =  <"oceanO 1 " , "ocean02 " , "ocean03" , 
"ocean04", "ocean05", "ocean06", 

"ocean07", "ocean08", "ocean09", 

"oceanlO", "wakeOla", "wakeOlb", 
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</boundary> 


ductOla" , "duct Olb" , "duct04a" , 
duct04b" , "duct05a" , "duct05b"> 


!  hole  surface  definitions 

<surface  name  =  "phantomOla  hole  boundary"> 
ijk_range  =  1,  1,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantomOla  hole  boundary"> 
ijk_range  =  3,  3,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantomOla  hole  boundary"> 
ijk_range  =  1,  3,  1,  1,  1,41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantomOla  hole  boundary"> 
ijk_range  =  1,  3,  31,  31,  1,  41, 

boundary^condition  —"cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantomOla  hole  boundary"> 
ijk_range  =  1,  3,  1,  31,  41,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom02a  hole  boundary"> 
ijk^range  =  1,  1,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom02a  hole  boundary"> 
ijk^range  =  3,  3,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom02a  hole  boundary"> 
ijk^range  =  1,  3,  1,  1,  1,41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom02a  hole  boundary"> 
ijk^range  =  1,  3,  31,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom02a  hole  boundary"> 
ijk^range  =  1,  3,  1,  31,  41,  41, 

boundary^condition  ="cut". 
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surface  normal  = 
</surf ace> 


<surface  name  =  "phantom03a  hole  boundary"> 
ijk^range  =  1,  1,  1,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom03a  hole  boundary"> 
ijk_range  =  3,  3,  1,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom03a  hole  boundary"> 
ijk_range  =  1,  3,  1,  1,  1,41, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom03a  hole  boundary"> 
ijk  range  =  1,  3,  31,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom03a  hole  boundary"> 
ijk  range  =  1,  3,  1,  31,  41,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom04a  hole  boundary"> 
ijk^range  =  1,  1,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom04a  hole  boundary"> 
ijk^range  =  3,  3,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom04a  hole  boundary"> 
ijk_range  =  1,  3,  1,  1,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom04a  hole  boundary"> 
ijk_range  =  1,  3,  31,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom04a  hole  boundary"> 
ijk_range  =  1,  3,  1,  31,  41,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 
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<surface  name  =  "phantomOlb  hole  boundary"> 
ijk_range  =  1,  1,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantomOlb  hole  boundary"> 
ijk_range  =  3,  3,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantomOlb  hole  boundary"> 
ijk_range  =  1,  3,  1,  1,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantomOlb  hole  boundary"> 
ijk^range  =  1,  3,  31,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantomOlb  hole  boundary"> 
ijk^range  =  1,  3,  1,  31,  41,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom02b  hole  boundary"> 
ijk^range  =  1,  1,  1,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom02b  hole  boundary"> 
ijk^range  =  3,  3,  1,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "tijk", 

</surf ace> 

<surface  name  =  "phantom02b  hole  boundary"> 
ijk^range  =  1,  3,  1,  1,  1,41, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom02b  hole  boundary"> 
ijk^range  =  1,  3,  31,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom02b  hole  boundary"> 
ijk^range  =  1,  3,  1,  31,  41,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom03b  hole  boundary"> 
ijk_range  =  1,  1,  1,  31,  1,  41, 
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boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom03b  hole  boundary"> 
ijk^range  =  3,  3,  1,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom03b  hole  boundary"> 
ijk^range  =  1,  3,  1,  1,  1,41, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom03b  hole  boundary"> 
ijk^range  =  1,  3,  31,  31,  1,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom03b  hole  boundary"> 
ijk_range  =  1,  3,  1,  31,  41,  41, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom04b  hole  boundary"> 
ijk_range  =  1,  1,  1,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom04b  hole  boundary"> 
ijk^range  =  3,  3,  1,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom04b  hole  boundary"> 
ijk_range  =  1,  3,  1,  1,  1,41, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom04b  hole  boundary"> 
ijk_range  =  1,  3,  31,  31,  1,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom04b  hole  boundary"> 
ijk_range  =  1,  3,  1,  31,  41,  41, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "phantom05  hole  boundary"> 
ijk_range  =  1,  1,  1,  3,  1,  2, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 
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<surface  name  =  "phantom05  hole  boundary"> 
ijk^range  =  2,  2,  1,  3,  1,  2, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom05  hole  boundary"> 
ijk^range  =  1,  2,  1,  1,  1,  2, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom05  hole  boundary"> 
ijk^range  =  1,  2,  3,  3,  1,  2, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "phantom05  hole  boundary"> 
ijk_range  =  1,  2,  1,  3,  1,  1, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "phantom05  hole  boundary"> 
ijk_range  =  1,  2,  1,  3,  2,  2, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "ductOla  hole  boundary"> 

ijk_range  =  7,  7,  1,  37,  1,  121, 

boundary^condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "ductOla  hole  boundary"> 

ijk_range  =  7,  21,  37,  37,  1,  121, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "duct04a  hole  boundary"> 

ijk_range  =  1,  53,  3,  3,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "duct05a  hole  boundary"> 

ijk_range  =  1,  15,  59,  59,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "wakeOla  hole  boundary"> 

ijk  range  =  1,  1,  1,  29,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk". 
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</surf ace> 

<surface  name  =  "wakeOla  hole  boundary"> 

ijk_range  =  36,  36,  1,  29,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "wakeOla  hole  boundary"> 

ijk_range  =  1,  36,  29,  29,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "ductOlb  hole  boundary"> 

ijk^range  =  7,  7,  1,  37,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "ductOlb  hole  boundary"> 

ijk^range  =  7,  21,  37,  37,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "duct04b  hole  boundary"> 

ijk^range  =  1,  53,  3,  3,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "duct05b  hole  boundary"> 

ijk^range  =  1,  15,  59,  59,  1,  121, 

boundary^condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "wakeOlb  hole  boundary"> 

ijk_range  =  1,  1,  1,  29,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "-ijk", 

</surf ace> 

<surface  name  =  "wakeOlb  hole  boundary"> 

ijk_range  =  36,  36,  1,  29,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 

<surface  name  =  "wakeOlb  hole  boundary"> 

ijk_range  =  1,  36,  29,  29,  1,  121, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 

</surf ace> 


<surface  name  =  "bargeOl  hole  boundary"> 
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1,  61, 


ijk_range  =  1,  107,  17,  17, 

boundary  condition  ="cut", 
surface  normal  =  "+ijk", 
</surf ace> 


!  outer  boundary  definition 

<boundary  name  =  "propellerO la  outer  boundary"> 
parent  grid  =  "propellerO la" , 

</boundary> 

<boundary  name  =  "propeller02a  outer  boundary"> 
parent  grid  =  "propeller02a" , 

</boundary> 

<boundary  name  =  "propeller03a  outer  boundary"> 
parent  grid  =  "propeller03a" , 

</boundary> 

<boundary  name  =  "propeller04a  outer  boundary"> 
parent  grid  =  "propeller04a" , 

</boundary> 

<boundary  name  =  "tipOla  outer  boundary"> 
parent  grid  =  "tipOla", 

</boundary> 

<boundary  name  =  "tip02a  outer  boundary"> 
parent  grid  =  "tip02a", 

</boundary> 

<boundary  name  =  "tip03a  outer  boundary"> 
parent_grid  =  "tip03a", 

</boundary> 

<boundary  name  =  "tip04a  outer  boundary"> 
parent  grid  =  "tip04a", 

</boundary> 

<boundary  name  =  "ductOla  outer  boundary"> 
parent  grid  =  "ductOla", 

</boundary> 

<boundary  name  =  "duct02a  outer  boundary"> 
parent  grid  =  "duct02a", 

</boundary> 

<boundary  name  =  "duct03a  outer  boundary"> 
parent_grid  =  "duct03a", 

</boundary> 

<boundary  name  =  "duct04a  outer  boundary"> 
parent  grid  =  "duct04a", 

</boundary> 

<boundary  name  =  "duct05a  outer  boundary"> 
parent_grid  =  "duct05a", 

</boundary> 

<boundary  name  =  "wakeOla  outer  boundary"> 
parent  grid  =  "wakeOla", 

</boundary> 

<boundary  name  =  "propellerO lb  outer  boundary"> 
parent  grid  =  "propellerOlb" , 

</boundary> 

<boundary  name  =  "propeller02b  outer  boundary"> 
parent  grid  =  "propeller02b" , 
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</boundary> 

<boundary  name  =  "propeller03b  outer  boundary"> 
parent  grid  =  "propeller03b" , 

</boundary> 

<boundary  name  =  "propeller04b  outer  boundary"> 
parent  grid  =  "propeller04b" , 

</boundary> 

<boundary  name  =  "tipOlb  outer  boundary"> 
parent_grid  =  "tipOlb", 

</boundary> 

<boundary  name  =  "tip02b  outer  boundary"> 
parent  grid  =  "tip02b", 

</boundary> 

<boundary  name  =  "tip03b  outer  boundary"> 
parent_grid  =  "tip03b", 

</boundary> 

<boundary  name  =  "tip04b  outer  boundary"> 
parent  grid  =  "tip04b", 

</boundary> 

<boundary  name  =  "ductOlb  outer  boundary"> 
parent  grid  =  "ductOlb", 

</boundary> 

<boundary  name  =  "duct02b  outer  boundary"> 
parent  grid  =  "duct02b", 

</boundary> 

<boundary  name  =  "duct03b  outer  boundary"> 
parent_grid  =  "duct03b", 

</boundary> 

<boundary  name  =  "duct04b  outer  boundary"> 
parent  grid  =  "duct04b", 

</boundary> 

<boundary  name  =  "duct05b  outer  boundary"> 
parent_grid  =  "duct05b", 

</boundary> 

<boundary  name  =  "wakeOlb  outer  boundary"> 
parent  grid  =  "wakeOlb", 

</boundary> 

<boundary  name  =  "oceanOl  outer  boundary"> 
parent  grid  =  "oceanOl", 

</boundary> 

<boundary  name  =  "ocean02  outer  boundary"> 
parent  grid  =  "ocean02", 

</boundary> 

<boundary  name  =  "ocean03  outer  boundary"> 
parent  grid  =  "ocean03", 

</boundary> 

<boundary  name  =  "ocean04  outer  boundary"> 
parent  grid  =  "ocean04", 

</boundary> 

<boundary  name  =  "ocean05  outer  boundary"> 
parent  grid  =  "ocean05", 

</boundary> 

<boundary  name  =  "ocean06  outer  boundary"> 
parent  grid  =  "ocean06", 

</boundary> 

<boundary  name  =  "ocean07  outer  boundary"> 
parent  grid  =  "ocean07". 
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</boundary> 

<boundary  name  =  "ocean08  outer 
parent  grid  =  "ocean08", 
</boundary> 

<boundary  name  =  "ocean09  outer 
parent  grid  =  "ocean09", 
</boundary> 

<boundary  name  =  "oceanlO  outer 
parent  grid  =  "oceanlO", 
</boundary> 

<boundary  name  =  "oceanll  outer 
parent  grid  =  "oceanll", 
</boundary> 

<boundary  name  =  "oceanl2  outer 
parent  grid  =  "oceanl2", 
</boundary> 

<boundary  name  =  "oceanl3  outer 
parent  grid  =  "oceanl3", 
</boundary> 

<boundary  name  =  "oceanl4  outer 
parent  grid  =  "oceanl4", 
</boundary> 

<boundary  name  =  "oceanl5  outer 
parent  grid  =  "oceanl5", 
</boundary> 

<boundary  name  =  "oceanl6  outer 
parent  grid  =  "oceanl6", 
</boundary> 

<boundary  name  =  "oceanl7  outer 
parent  grid  =  "oceanl7", 
</boundary> 

<boundary  name  =  "oceanl8  outer 
parent  grid  =  "oceanl8", 
</boundary> 

<boundary  name  =  "barged  outer 
parent  grid  =  "bargeOl", 
</boundary> 


boundary"> 

boundary"> 

boundary"> 

boundary"> 

boundary"> 

boundary"> 

boundary"> 

boundary"> 

boundary"> 

boundary"> 

boundary"> 

boundary"> 


!  outer  boundary  surface  definition 


<surface  name  =  "propellerO la  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propellerO la" 
donor_i j k_range  =  61,  61,  1,  35,  1,  42, 

</surf ace> 

<surface  name  =  "propellerO la  outer  boundary"> 
ijk_range  =  62,  62,1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propellerOla" 
donor_ijk  range  =  2,2,1,  35,1,  42, 

</surf ace> 

<surface  name  =  "propellerOla  outer  boundary"> 
ijk_range  =  1,  62,  35,  35,  2,  41, 

</surf ace> 

<surface  name  =  "propellerOla  outer  boundary"> 
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1,  62, 


2,  35,  42,  42, 


ijk  range  = 
</surf ace> 


<surface  name  =  "propeller02a  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller02a" 
donor_i j k_range  =  61,  61,  1,  35,  1,  42, 

</surf ace> 

<surface  name  =  "propeller02a  outer  boundary"> 
ijk_range  =  62,  62,1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller02a" 
donor_ijk  range  =  2,2,1,  35,1,  42, 

</surf ace> 

<surface  name  =  "propeller02a  outer  boundary"> 
ijk_range  =  1,  62,  35,  35,  2,  41, 

</surf ace> 

<surface  name  =  "propeller02a  outer  boundary"> 
ijk_range  =  1,  62,  2,  35,  42,  42, 

</surf ace> 

<surface  name  =  "propeller03a  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller03a" 
donor_i j k_range  =  61,  61,  1,  35,  1,  42, 

</surf ace> 

<surface  name  =  "propeller03a  outer  boundary"> 
ijk_range  =  62,  62,1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller03a" 
donor_ijk  range  =  2,2,1,  35,1,  42, 

</surf ace> 

<surface  name  =  "propeller03a  outer  boundary"> 
ijk_range  =  1,  62,  35,  35,  2,  41, 

</surf ace> 

<surface  name  =  "propeller03a  outer  boundary"> 
ijk_range  =  1,  62,  2,  35,  42,  42, 

</surf ace> 

<surface  name  =  "propeller04a  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller04a" 
donor_i j k_range  =  61,  61,  1,  35,  1,  42, 

</surf ace> 

<surface  name  =  "propeller04a  outer  boundary"> 
ijk_range  =  62,  62,1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller04a" 
donor_ijk  range  =  2,2,1,  35,1,  42, 

</surf ace> 

<surface  name  =  "propeller04a  outer  boundary"> 
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ijk_range  =  1,  62,  35,  35,  2,  41, 

</surf ace> 

<surface  name  =  "propeller04a  outer  boundary"> 
ijk_range  =  1,  62,  2,  35,  42,  42, 

</surf ace> 


<surface  name  =  "tipOla  outer  boundary"> 
ijk_range  =  1,  1,  1,  4,  2, 5, 

</surf ace> 

<surface  name  =  "tipOla  outer  boundary"> 
ijk  range  =  29,  29,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tipOla  outer  boundary"> 
ijk  range  =  1,  29,  1,  1,  1,  5, 

</surf ace> 

<surface  name  =  "tipOla  outer  boundary"> 
ijk  range  =  1,  29,  4,  4,  1,  5, 

</surf ace> 

<surface  name  =  "tipOla  outer  boundary"> 
ijk  range  =  1,  29,  1,  4,  5,  5, 

</surf ace> 


<surface  name  =  "tip02a  outer  boundary"> 
ijk_range  =  1,  1,  1,  4,  2, 5, 

</surf ace> 

<surface  name  =  "tip02a  outer  boundary"> 
ijk  range  =  29,  29,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tip02a  outer  boundary"> 
ijk  range  =  1,  29,  1,  1,  1,  5, 

</surf ace> 

<surface  name  =  "tip02a  outer  boundary"> 
ijk  range  =  1,  29,  4,  4,  1,  5, 

</surf ace> 

<surface  name  =  "tip02a  outer  boundary"> 
ijk  range  =  1,  29,  1,  4,  5,  5, 

</surf ace> 


<surface  name  =  "tip03a  outer  boundary"> 
ijk_range  =  1,  1,  1,  4,  2, 5, 

</surf ace> 

<surface  name  =  "tip03a  outer  boundary"> 
ijk  range  =  29,  29,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tip03a  outer  boundary"> 
ijk_range  =  1,  29,  1,  1,  1,  5, 

</surf ace> 

<surface  name  =  "tip03a  outer  boundary"> 
ijk  range  =  1,  29,  4,  4,  1,  5, 

</surf ace> 

<surface  name  =  "tip03a  outer  boundary"> 
ijk  range  =  1,  29,  1,  4,  5,  5, 

</surf ace> 
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<surface  name  =  "tip04a  outer  boundary"> 
ijk_range  =  1,  1,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tip04a  outer  boundary"> 
ijk  range  =  29,  29,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tip04a  outer  boundary"> 
ijk  range  =  1,  29,  1,  1,  1,  5, 

</surf ace> 

<surface  name  =  "tip04a  outer  boundary"> 
ijk  range  =  1,  29,  4,  4,  1,  5, 

</surf ace> 

<surface  name  =  "tip04a  outer  boundary"> 
ijk  range  =  1,  29,  1,  4,  5,  5, 

</surf ace> 


<surface  name  =  "ductOla  outer  boundary"> 
ijk  range  =  1,  1,  1,  57,  1,  122, 

</surf ace> 

<surface  name  =  "ductOla  outer  boundary"> 
ijk  range  =  21,  21,  1,  34,  1,  122, 

</surf ace> 

<surface  name  =  "ductOla  outer  boundary"> 
ijk  range  =  21,  21,  36,  57,  1,  122, 

</surf ace> 

<surface  name  =  "ductOla  outer  boundary"> 
ijk  range  =  8,  21,  1,  1,  1,  122, 

</surf ace> 

<surface  name  =  "ductOla  outer  boundary"> 
ijk  range  =  1,  21,  57,  57,  1,  122, 

</surf ace> 

<surface  name  =  "ductOla  outer  boundary"> 
ijk_range  =  1,  21,  1,  57,  1,  1, 

boundary  condition  =  "periodic", 
donor^grid  =  "ductOla" 

donor_ijk  range  =  1,  21,  1,  57,  121,  121, 
</surf ace> 

<surface  name  =  "ductOla  outer  boundary"> 

ijk  range  =  1,  21,  1,  57,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "ductOla" 

donor_ijk  range  =  1,  21,  1,  57,  2,  2, 

</surf ace> 


<surface  name  =  "duct02a  outer  boundary"> 
ijk_range  =  66,  66,  1,  24,  1,  122, 

</surf ace> 

<surface  name  =  "duct02a  outer  boundary"> 
ijk  range  =  1,  66,  24,  24,  1,  122, 

</surf ace> 

<surface  name  =  "duct02a  outer  boundary"> 
ijk_range  =  1,  66,  1,  24,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct02a" 
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1,  66 


1,  24,  121,  121 


donor_ijk  range  = 

</surf ace> 

<surface  name  =  "duct02a  outer  boundary"> 

ijk  range  =  1,  66,  1,  24,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct02a" 

donor_ijk  range  =  1,  66,  1,  24,  2,  2, 

</surf ace> 


<surface  name  =  "duct03a  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  122, 

</surf ace> 

<surface  name  =  "duct03a  outer  boundary"> 
ijk_range  =53,  53,  1,  35,  1,  122, 

</surf ace> 

<surface  name  =  "duct03a  outer  boundary"> 
ijk_range  =  1,  53,  1,  1,  1,  122, 

</surf ace> 

<surface  name  =  "duct03a  outer  boundary"> 
ijk_range  =  1,  53,  1,  35,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct03a" 

donor_i j k^range  =  1,  53,  1,  35,  121,  121, 

</surf ace> 

<surface  name  =  "duct03a  outer  boundary"> 

ijk_range  =  1,  53,  1,  35,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct03a" 

donor_i j k_range  =  1,  53,  1,  35,  2,  2, 

</surf ace> 


<surface  name  =  "duct04a  outer  boundary"> 
ijk  range  =  1,  1,  1,  23,  1,  122, 

</surf ace> 

<surface  name  =  "duct04a  outer  boundary"> 
ijk_range  =53,  53,  1,  23,  1,  122, 

</surf ace> 

<surface  name  =  "duct04a  outer  boundary"> 
ijk_range  =  1,  53,  23,  23,  1,  122, 

</surf ace> 

<surface  name  =  "duct04a  outer  boundary"> 
ijk_range  =  1,  53,  1,  23,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct04a" 

donor_ijk  range  =  1,  53,  1,  23,  121,  121, 

</surf ace> 

<surface  name  =  "duct04a  outer  boundary"> 

ijk  range  =  1,  53,  1,  23,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct04a" 

donor_i j k_range  =  1,  53,  1,  23,  2,  2, 

</surf ace> 


<surface  name  =  "duct05a  outer  boundary"> 
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1, 


1, 


1,  56, 


1,  122, 


ijk  range  = 

</surf ace> 

<surface  name  =  "duct05a  outer  boundary"> 
ijk_range  =  1,  1,  58,  79,  1,  122, 

</surf ace> 

<surface  name  =  "duct05a  outer  boundary"> 
ijk  range  =  21,  21,  1,  79,  1,  122, 

</surf ace> 

<surface  name  =  "duct05a  outer  boundary"> 
ijk  range  =  1,  21,  79,  79,  1,  122, 

</surf ace> 

<surface  name  =  "duct05a  outer  boundary"> 
ijk_range  =  1,  21,  1,  79,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct05a" 

donor_ijk  range  =  1,  21,  1,  79,  121,  121, 
</surf ace> 

<surface  name  =  "duct05a  outer  boundary"> 

ijk  range  =  1,  21,  1,  79,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct05a" 

donor_ijk  range  =  1,  21,  1,  79,  2,  2, 

</surf ace> 


<surface  name  =  "wakeOla  outer  boundary"> 
ijk  range  =  1,  1,  1,  41,  1,  122, 

</surf ace> 

<surface  name  =  "wakeOla  outer  boundary"> 
ijk  range  =  40,  40,  1,  41,  1,  122, 

</surf ace> 

<surface  name  =  "wakeOla  outer  boundary"> 
ijk  range  =  1,  40,  41,  41,  1,  122, 

</surf ace> 

<surface  name  =  "wakeOla  outer  boundary"> 
ijk_range  =  1,  40,  1,  41,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "wakeOla" 

donor_ijk  range  =  1,  40,  1,  41,  121,  121, 

</surf ace> 

<surface  name  =  "wakeOla  outer  boundary"> 

ijk  range  =  1,  40,  1,  41,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "wakeOla" 

donor_ijk  range  =  1,  40,  1,  41,  2,  2, 

</surf ace> 


<surface  name  =  "propellerO lb  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propellerO lb" 
donor_i j k_range  =  61,  61,  1,  35,  1,  42, 

</surf ace> 

<surface  name  =  "propellerO lb  outer  boundary"> 
ijk_range  =  62,  62,  1,  35,  1,  42, 

boundary  condition  =  "periodic". 
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donor_grid  =  "propellerOlb" 
donor_ijk  range  =  2,  2,  1,  35,  1,  42, 

</surf ace> 

<surface  name  =  "propellerOlb  outer  boundary"> 
ijk_range  =  1,  62,  35,  35,  2,  41, 

</surf ace> 

<surface  name  =  "propellerOlb  outer  boundary"> 
ijk_range  =  1,  62,  2,  35,  42,  42, 

</surf ace> 

<surface  name  =  "propeller02b  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller02b" 
donor_i j k_range  =  61,  61,  1,  35,  1,  42, 

</surf ace> 

<surface  name  =  "propeller02b  outer  boundary"> 
ijk_range  =  62,  62,1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller02b" 
donor_ijk  range  =  2,2,1,  35,1,  42, 

</surf ace> 

<surface  name  =  "propeller02b  outer  boundary"> 
ijk_range  =  1,  62,  35,  35,  2,  41, 

</surf ace> 

<surface  name  =  "propeller02b  outer  boundary"> 
ijk_range  =  1,  62,  2,  35,  42,  42, 

</surf ace> 

<surface  name  =  "propeller03b  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller03b" 
donor_i j k_range  =  61,  61,  1,  35,  1,  42, 

</surf ace> 

<surface  name  =  "propeller03b  outer  boundary"> 
ijk_range  =  62,  62,1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller03b" 
donor_ijk  range  =  2,2,1,  35,1,  42, 

</surf ace> 

<surface  name  =  "propeller03b  outer  boundary"> 
ijk_range  =  1,  62,  35,  35,  2,  41, 

</surf ace> 

<surface  name  =  "propeller03b  outer  boundary"> 
ijk_range  =  1,  62,  2,  35,  42,  42, 

</surf ace> 


<surface  name  =  "propeller04b  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller04b" 
donor_i j k_range  =  61,  61,  1,  35,  1,  42, 

</surf ace> 
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<surface  name  =  "propeller04b  outer  boundary"> 
ijk_range  =  62,  62,1,  35,  1,  42, 

boundary  condition  =  "periodic", 
donor_grid  =  "propeller04b" 
donor_ijk  range  =  2,2,1,  35,1,  42, 

</surf ace> 

<surface  name  =  "propeller04b  outer  boundary"> 
ijk_range  =  1,  62,  35,  35,  2,  41, 

</surf ace> 

<surface  name  =  "propeller04b  outer  boundary"> 
ijk_range  =  1,  62,  2,  35,  42,  42, 

</surf ace> 


<surface  name  =  "tipOlb  outer  boundary"> 
ijk_range  =  1,  1,  1,  4,  2, 5, 

</surf ace> 

<surface  name  =  "tipOlb  outer  boundary"> 
ijk  range  =  29,  29,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tipOlb  outer  boundary"> 
ijk_range  =  1,  29,  1,  1,  1,  5, 

</surf ace> 

<surface  name  =  "tipOlb  outer  boundary"> 
ijk  range  =  1,  29,  4,  4,  1,  5, 

</surf ace> 

<surface  name  =  "tipOlb  outer  boundary"> 
ijk  range  =  1,  29,  1,  4,  5,  5, 

</surf ace> 


<surface  name  =  "tip02b  outer  boundary"> 
ijk_range  =  1,  1,  1,  4,  2, 5, 

</surf ace> 

<surface  name  =  "tip02b  outer  boundary"> 
ijk  range  =  29,  29,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tip02b  outer  boundary"> 
ijk  range  =  1,  29,  1,  1,  1,  5, 

</surf ace> 

<surface  name  =  "tip02b  outer  boundary"> 
ijk  range  =  1,  29,  4,  4,  1,  5, 

</surf ace> 

<surface  name  =  "tip02b  outer  boundary"> 
ijk  range  =  1,  29,  1,  4,  5,  5, 

</surf ace> 


<surface  name  =  "tip03b  outer  boundary"> 
ijk_range  =  1,  1,  1,  4,  2, 5, 

</surf ace> 

<surface  name  =  "tip03b  outer  boundary"> 
ijk  range  =  29,  29,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tip03b  outer  boundary"> 
ijk_range  =  1,  29,  1,  1,  1,  5, 

</surf ace> 
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<surface  name  =  "tip03b  outer  boundary"> 
ijk  range  =  1,  29,  4,  4,  1,  5, 

</surf ace> 

<surface  name  =  "tip03b  outer  boundary"> 
ijk  range  =  1,  29,  1,  4,  5,  5, 

</surf ace> 


<surface  name  =  "tip04b  outer  boundary"> 
ijk_range  =  1,  1,  1,  4,  2, 5, 

</surf ace> 

<surface  name  =  "tip04b  outer  boundary"> 
ijk  range  =  29,  29,  1,  4,  2,  5, 

</surf ace> 

<surface  name  =  "tip04b  outer  boundary"> 
ijk  range  =  1,  29,  1,  1,  1,  5, 

</surf ace> 

<surface  name  =  "tip04b  outer  boundary"> 
ijk  range  =  1,  29,  4,  4,  1,  5, 

</surf ace> 

<surface  name  =  "tip04b  outer  boundary"> 
ijk  range  =  1,  29,  1,  4,  5,  5, 

</surf ace> 


<surface  name  =  "ductOlb  outer  boundary"> 
ijk  range  =  1,  1,  1,  57,  1,  122, 

</surf ace> 

<surface  name  =  "ductOlb  outer  boundary"> 
ijk  range  =  21,  21,  1,  34,  1,  122, 

</surf ace> 

<surface  name  =  "ductOlb  outer  boundary"> 
ijk  range  =  21,  21,  36,  57,  1,  122, 

</surf ace> 

<surface  name  =  "ductOlb  outer  boundary"> 
ijk  range  =  8,  21,  1,  1,  1,  122, 

</surf ace> 

<surface  name  =  "ductOlb  outer  boundary"> 
ijk  range  =  1,  21,  57,  57,  1,  122, 

</surf ace> 

<surface  name  =  "ductOlb  outer  boundary"> 
ijk_range  =  1,  21,  1,  57,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "ductOlb" 

donor_ijk  range  =  1,  21,  1,  57,  121,  121, 

</surf ace> 

<surface  name  =  "ductOlb  outer  boundary"> 

ijk  range  =  1,  21,  1,  57,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "ductOlb" 

donor_ijk  range  =  1,  21,  1,  57,  2,  2, 

</surf ace> 


<surface  name  =  "duct02b  outer  boundary"> 
ijk_range  =66,  66,  1,  24,  1,  122, 

</surf ace> 
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<surface  name  =  "duct02b  outer  boundary"> 
ijk  range  =  1,  66,  24,  24,  1,  122, 

</surf ace> 

<surface  name  =  "duct02b  outer  boundary"> 
ijk_range  =  1,  66,  1,  24,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct02b" 

donor_ijk  range  =  1,  66,  1,  24,  121,  121, 

</surf ace> 

<surface  name  =  "duct02b  outer  boundary"> 

ijk  range  =  1,  66,  1,  24,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct02b" 

donor_ijk  range  =  1,  66,  1,  24,  2,  2, 

</surf ace> 


<surface  name  =  "duct03b  outer  boundary"> 
ijk_range  =  1,  1,  1,  35,  1,  122, 

</surf ace> 

<surface  name  =  "duct03b  outer  boundary"> 
ijk_range  =53,  53,  1,  35,  1,  122, 

</surf ace> 

<surface  name  =  "duct03b  outer  boundary"> 
ijk_range  =  1,  53,  1,  1,  1,  122, 

</surf ace> 

<surface  name  =  "duct03b  outer  boundary"> 
ijk_range  =  1,  53,  1,  35,  1,  1, 

boundary  condition  =  "periodic", 
donor_jgrid  =  "duct03b" 

donor_ijk  range  =  1,  53,  1,  35,  121,  121, 
</surf ace> 

<surface  name  =  "duct03b  outer  boundary"> 

ijk_range  =  1,  53,  1,  35,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct03b" 

donor_i j k_range  =  1,  53,  1,  35,  2,  2, 
</surf ace> 


<surface  name  =  "duct04b  outer  boundary"> 
ijk  range  =  1,  1,  1,  23,  1,  122, 

</surf ace> 

<surface  name  =  "duct04b  outer  boundary"> 
ijk_range  =53,  53,  1,  23,  1,  122, 

</surf ace> 

<surface  name  =  "duct04b  outer  boundary"> 
ijk_range  =  1,  53,  23,  23,  1,  122, 

</surf ace> 

<surface  name  =  "duct04b  outer  boundary"> 
ijk_range  =  1,  53,  1,  23,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct04b" 

donor_i j k^range  =  1,  53,  1,  23,  121,  121, 

</surf ace> 

<surface  name  =  "duct04b  outer  boundary"> 

ijk  range  =  1,  53,  1,  23,  122,  122, 
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boundary  condition  =  "periodic", 
donor_grid  =  "duct04b" 

donor^i j k_range  =  1,  53,  1,  23,  2,  2, 

</surf ace> 


<surface  name  =  "duct05b  outer  boundary"> 
ijk_range  =  1,  1,  1,  56,  1,  122, 

</surf ace> 

<surface  name  =  "duct05b  outer  boundary"> 
ijk_range  =  1,  1,  58,  79,  1,  122, 

</surf ace> 

<surface  name  =  "duct05b  outer  boundary"> 
ijk  range  =  21,  21,  1,  79,  1,  122, 

</surf ace> 

<surface  name  =  "duct05b  outer  boundary"> 
ijk  range  =  1,  21,  79,  79,  1,  122, 

</surf ace> 

<surface  name  =  "duct05b  outer  boundary"> 
ijk_range  =  1,  21,  1,  79,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct05b" 

donor_ijk  range  =  1,  21,  1,  79,  121,  121, 

</surf ace> 

<surface  name  =  "duct05b  outer  boundary"> 

ijk  range  =  1,  21,  1,  79,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "duct05b" 

donor_ijk  range  =  1,  21,  1,  79,  2,  2, 

</surf ace> 


<surface  name  =  "wakeOlb  outer  boundary"> 
ijk  range  =  1,  1,  1,  41,  1,  122, 

</surf ace> 

<surface  name  =  "wakeOlb  outer  boundary"> 
ijk  range  =  40,  40,  1,  41,  1,  122, 

</surf ace> 

<surface  name  =  "wakeOlb  outer  boundary"> 
ijk  range  =  1,  40,  41,  41,  1,  122, 

</surf ace> 

<surface  name  =  "wakeOlb  outer  boundary"> 
ijk_range  =  1,  40,  1,  41,  1,  1, 

boundary  condition  =  "periodic", 
donor_grid  =  "wakeOlb" 

donor_ijk  range  =  1,  40,  1,  41,  121,  121, 

</surf ace> 

<surface  name  =  "wakeOlb  outer  boundary"> 

ijk  range  =  1,  40,  1,  41,  122,  122, 

boundary  condition  =  "periodic", 
donor_grid  =  "wakeOlb" 

donor__ijk  range  =  1,  40,  1,  41,  2,  2, 

</surf ace> 


<surface  name  =  "oceanOl  outer  boundary"> 
ijk_range  =  20,  20,  1,  116,  1,  95, 


C-25 


</surf ace> 

<surface  name  =  "oceanOl  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 


"ocean02  outer  boundary"> 

=  1,  1,  1,  H6,  1,  95, 

"ocean02  outer  boundary"> 

=  20,  20,  1,  116,  1,  95, 

"ocean02  outer  boundary"> 

=  1,  20,  1,  1,  1,  95, 


<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 


"ocean03  outer  boundary"> 

=  1,  1,  1,  H6,  1,  95, 

"ocean03  outer  boundary"> 

=  20,  20,  1,  116,  1,  95, 

"ocean03  outer  boundary"> 

=  1,  20,  1,  1,  1,  95, 


<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 


"ocean04  outer  boundary"> 

=  1,  1,  1,  H6,  1,  95, 

"ocean04  outer  boundary"> 

=  20,  20,  1,  116,  1,  95, 

"ocean04  outer  boundary"> 

=  1,  20,  1,  1,  1,  95, 


<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 


"ocean05  outer  boundary"> 

=  1,  1,  1,  H6,  1,  95, 

"ocean05  outer  boundary"> 

=  20,  20,  1,  116,  1,  95, 

"ocean05  outer  boundary"> 

=  1,  20,  1,  1,  1,  95, 


<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 
ijk  range 
</surf ace> 
<surface  name  = 


"ocean06  outer  boundary"> 

=  1,  1,  1,  H6,  1,  95, 

"ocean06  outer  boundary"> 

=  20,  20,  1,  116,  1,  95, 

"ocean06  outer  boundary"> 
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ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "ocean07  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "ocean07  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "ocean07  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "ocean08  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "ocean08  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "ocean08  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "ocean09  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "ocean09  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "ocean09  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "oceanlO  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanlO  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanlO  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "oceanll  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanll  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanll  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 
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<surface  name  =  "oceanl2  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl2  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl2  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "oceanl3  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl3  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl3  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "oceanl4  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl4  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl4  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "oceanl5  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl5  outer  boundary"> 
ijk_range  =20,  20,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl5  outer  boundary"> 
ijk_range  =  1,  20,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "oceanl6  outer  boundary"> 
ijk_range  =  1,  1,  1,  116,  1,  95, 

</surf ace> 

<surface  name  =  "oceanl6  outer  boundary"> 
ijk_range  =  1,  16,  1,  1,  1,  95, 

</surf ace> 


<surface  name  =  "oceanl7  outer  boundary"> 
ijk_range  =  1,  151,  37,  37,  1,  33, 

</surf ace> 
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<surface  name  =  "oceanl7  outer  boundary"> 
ijk_range  =  1,  151,  1,  37,  33,  33, 

</surf ace> 


<surface  name  =  "oceanl8  outer  boundary"> 
ijk_range  =  1,  151,  37,  37,  1,  33, 

</surf ace> 

<surface  name  =  "oceanl8  outer  boundary"> 
ijk_range  =  1,  151,  1,  37,  1,  1, 

</surf ace> 


<surface  name  =  "barged  outer  boundary"> 
ijk_range  =  1,  107,  34,  34,  1,  61, 

</surf ace> 
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